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Major  Department:  Chemistry 

The  infrared  spectra  of  uracil  and  methylated  uracils  have  been 
interpreted  by  carrying  out  theoretical  predictions  of  these  spectra 
for  such  biomolecular  systems.  The  infrared  spectra  for  these  mole- 
cules isolated  in  noble  gas  matrices  are  interpreted  by  comparison 
with  the  predicted  spectra.  These  predictions  are  made  using  force 
constants  and  intensity  parameters  (Atomic  Polar  Tensors)  from 
quantum  mechanical  ab  initio  calculations  with  a 4— 31G  basis  set 
for  the  in-plane  modes  of  these  uracils. 

The  infrared  frequencies  and  intensities  of  the  methylated  uracils 
are  predicted  transferring  vibrational  parameters  from  other  molecules. 
The  absolute  intensities  predicted  for  these  molecules  are  compared 
with  the  experimental  relative  intensities,  which  were  obtained  in 
matrices.  Comparisons  of  ab  initio  and  semi-empirical  force  constants 
are  discussed,  as  well  as  intensity  parameters. 


xi 


The  comparison  of  the  transferred  atomic  polar  tensors  and  of  the 
force  constants  from  quantum  mechanical  calculations  is  discussed  for 
each  of  the  following  molecules:  1-methyluracil,  l-methyl-3-D-uracil , 
3-methyluracil,  and  1,3-dimethyluracil.  This  comparison  shows  that 
atomic  polar  tensors  transferred  from  uracil,  acetone  and  N-methyl- 
acetamide  predict  the  observed  intensities  of  the  methylated  uracils 
within  the  usual  "factor-of-two" . Also,  the  scaled  ab  initio  force 
constants  reproduce  the  experimental  wavenumbers  of  a molecule 
within  about  20  cm  ^ . 
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CHAPTER  1 
INTRODUCTION 

In  recent  years,  vibrational  spectroscopy  has  proven  to  be  a 
powerful  tool  in  biophysical  research  on  nucleic  acids  (1-6).  In  par- 
ticular, infrared  and  Raman  spectra  have  played  a critical  role  in  the 
study  of  molecular  structures  and  dynamical  properties  of  nucleosides, 
nucleotides,  and  the  constituent  bases  of  nucleic  acids  (2,4,6). 

However,  the  full  power  of  infrared  and  Raman  spectroscopy  depends  not 
only  on  obtaining  spectra  of  such  biomolecules,  but  on  the  most  accurate 
interpretation  of  their  spectra.  In  order  to  properly  accomplish  this 
purpose,  it  is  essential  first  to  obtain,  and  then  to  understand,  the 
details  of  the  spectra  for  the  isolated  non-interacting  monomeric 
molecules.  It  is  then  required  that  these  experimental  data  be  comple- 
mented with  theoretical  calculations.  This  procedure  leads  to  the 
establishment  of  correlations  between  spectral  changes  which  occur  from 
the  isolated  molecule  to  molecules  placed  in  more  and  more  strongly 
interacting  environments. 

It  is  well  known  that  nucleic  acids  are  important  constituents  of 
all  living  cells  because  they  are  the  biomolecules  which  carry  and 
transmit  genetic  information  used  in  the  synthesis  of  the  cellular 
proteins  (7,8).  Depending  upon  the  molecular  structure  of  their  con- 
stituent sugars,  nucleic  acids  may  be  divided  into  two  distinct  classes, 
ribonucleic  acid  (RNA)  or  deoxyribonucleic  acid  (DNA) . In  RNA  the  sugar 

is  D-ribofuranose  and  in  DNA  is  2-deoxy-D-ribofuranose.  Both  RNA  and 
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DNA  are  the  result  of  successive  linking  of  alternating  simple  units 
consisting  of  phosphoric  acid,  sugar  residues,  and  purine  or  pyrimidine 
bases  (7).  These  simple  units,  which  are  called  nucleotides,  have  the 
general  structure  given  in  Fig.  1-1 (a)  for  a ribonucleotide.  For  a 
deoxyribonucleotide  the  OH  group  attached  to  C2'  in  the  ribose  ring  is 
replaced  by  an  H atom.  Thus,  a nucleotide  is  a phosphate  ester  of  a 
nucleoside,  like  uridine  (Fig.  1-1 (b)). 

It  is  known  that  the  genetic  code  consists  of  64  triplet  combina- 
tions (codons)  of  nucleotides,  as  shown  in  Fig.  1-2.  Of  the  64 
possible  combinations,  61  of  them  specify  particular  amino  acids  (7). 
The  sequence  of  codons  in  messenger-RNA  (mRNA)  is  translated  to  the 
sequence  of  amino  acids  in  the  protein  biosynthesis,  making  possible 
the  biological  process  of  transmission  of  genetic  information.  If  one 
of  the  bases  in  a codon  is  substituted  by  another,  (e.g.,  CGU  is 
changed  to  CAG  (Fig.  1-2)),  then  a mutation  occurs  because  this  new 
codon  would  code  for  a different  amino  acid  in  the  cell  (8) . 

As  seen  in  Fig.  1-2  one  very  important  nucleoside  in  the  genetic 
code  is  uridine  (base  = uracil)  which  is  present  in  37  of  the  64 
possible  codons.  Because  of  the  central  role  that  uridine  plays  on 
the  transmission  of  genetic  information,  it  is  essential  to  understand 
its  molecular  structure  in  complete  detail.  For  this  it  is  important 
to  carry  out  studies  of  its  isolated  molecular  components  in  order  to 
understand  their  interactions  with  other  molecules.  Further,  the 
systematic  study  of  the  structural  and  dynamical  properties  of  these 
molecules  represents  an  interesting  problem  in  its  own  rights. 

We  shall  be  concerned  in  this  dissertation  with  the  study  of  the 
uracil  base  and  some  of  its  methylated  derivatives.  Specifically  we 
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Figure  1—1.  Molecular  structures  of  nucleic  acid  constituents. 

(a)  General  structure  of  a ribonucleotide.  When 
the  OH  group  in  position  C2’  of  the  ring  system 

is  replaced  by  H,  then  it  is  a deoxyribonucleotide. 

(b)  Structure  of  uridine  (a  ribonucleoside) . 

(c)  Usual  pyrimidine  (uracil  and  cytosine)  and 
purine  (adenine  and  guanine)  bases.  To  form  the 
nucleotide  or  nucleoside,  the  N1  atom  of  the 
pyrimidine  or  the  N9  atom  of  the  purine  is 
attached  to  the  Cl'  atom  of  the  ribose  sugar. 
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shall  be  concerned  with  the  infrared  frequencies  and  intensities  of 
isolated  uracils,  and  in  particular  we  shall  carry  out  normal  coordinate 
analysis  calculations  for  the  in-plane  vibrations  of  these  molecules, 
using  these  together  with  the  dipole  moment  derivatives  to  predict 
quantitatively  their  infrared  spectra.  Then  the  calculated  and  experi- 
mental spectra  for  the  isolated  molecules  will  be  compared  so  that  more 
reliable  assignments  can  be  made  for  the  observed  spectra. 

We  believe  the  study  of  the  methylated  uracils  to  be  of  consider- 
able importance;  in  particular  1-methyluracil,  because  this  molecule  is 
more  representative  spectroscopically  of  the  base  in  the  uridylate 
nucleotide  or  uridine  nucleoside  than  uracil  itself,  since  the  heavy 
methyl  group  has  almost  the  same  effect  on  the  uracil  base  vibrations 
as  does  the  ribose  sugar.  The  study  of  the  other  methylated  uracils  is 
also  of  importance  because  they  provide  useful  tests  of  the  procedure 
used  for  assignment  of  1-methyluracil  and  their  spectra  provide  insight 
as  to  changes  expected  when  the  structure  changes  slightly.  Also,  the 
changes  in  the  spectrum  of  1 ,3-dimethyluracil  from  matrix-isolated 
molecule  to  solid  provide  some  insight  as  to  the  effect  of  stacking 
interactions,  since  this  molecule  is  free  of  hydrogen  bonding,  which 
dominates  the  same  changes  in  the  spectra  of  the  other  molecules. 

Historical  Perspective 

A considerable  number  of  studies  dealing  with  the  infrared  and 
Raman  spectra  of  uracil,  its  isotopic  derivatives  (deuterated,  and 

■*-®0)  , and  its  methylated  derivatives  have  been  reported  (1-6,8-20). 

We  shall  review  here  briefly  some  of  these  studies,  in  particular  those 
in  which  an  interpretation  of  the  spectra  has  been  made  by  means  of 
normal  coordinate  analyses  (NCAs) . 
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An  extensive  study  by  Angell  (9)  reported  the  infrared  spectra  in 
the  solid  state  of  uracil,  cytosine,  adenine,  and  guanine  bases,  of 
ribose  and  deoxyribose  nucleosides  derived  from  these  bases,  and  of  some 
mononucleotides.  He  found  that  in  the  case  of  uracil  and  cytosine,  only 
small  changes  occurred  in  the  part  of  the  spectrum  from  the  base  residue 
when  a sugar  group  was  added.  He  observed  a number  of  strong  bands  in 
the  spectrum  of  the  nucleoside  between  900  - 1200  cm  ^ which  are 
characteristic  of  the  sugar  residue.  Upon  addition  of  phosphate  group 
to  uridine  there  was  relatively  little  change  in  the  bands  due  to  the 
uridine.  He  also  found  that  the  spectra  of  uridine,  cytidine,  adeno- 
sine and  thymidine  do  not  change  from  that  for  the  solid  when  the 
compounds  are  dissolved  in  water.  According  to  Angell  (9)  this  finding 
is  important  when  one  attempts  to  extend  the  results  from  simple 
compounds  in  the  solid  state  to  the  spectra  of  nucleic  acids  that 
usually  exist  in  a highly  hydrated  state. 

The  rather  exhaustive  Raman  spectral  studies  of  Lord  and  Thomas 
(10,11)  showed  that  uracil  and  its  derivatives  in  any  given  environment 
have  similar  Raman  lines  in  the  1500  - 1800  cm  1 region  regardless  of 
the  substituents  at  position  (see  Fig.  l-l(c)).  They  interpreted 
the  Raman  spectra  of  these  molecules  in  terms  of  changes  in  molecular 
structure  produced  by  changes  in  pH,  pD  and  state  of  aggregation. 

In  this  study  an  analogy  was  established  between  the  observed  Raman 
spectra  of  uridine,  1-methyluracil  and  1,3-dimethyluracil.  They  found 
that  the  Raman  spectra  as  a function  of  pH  of  these  methylated  uracils 
resemble  those  of  uridine. 

In  the  study  by  Susi  and  Ard  (12)  laser-Raman  and  infrared  spectra 
of  polycrystalline  uracil,  N,N-dideutero-,  C,C-dideutero-,  and  perdeu- 
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tero  uracil  were  obtained.  In  order  to  make  assignments  of  these 
spectra,  they  carried  out  NCAs  for  the  planar  modes,  which  were  des- 
cribed in  terms  of  potential  energy  distributions  and  Cartesian  dis- 
placements. For  this  calculation  they  constructed  a 23  parameter 
valence  force  field  by  a least-squares  refinement  procedure,  starting 
by  transferring  force  constants  from  other  molecules,  for  example  ace- 
tone, to  uracil.  This  force  field  reproduced  83  observed  frequencies 
for  uracil  and  its  deuterated  derivatives,  with  an  average  error  of 
8 cm-  . 

More  NCAs  were  carried  out  for  pyrimidine  derivatives  by  Susi  and 
Ard  (13).  These  included  cytosine,  thymine,  1-methyluracil , and  some 
deutero  derivatives  of  these  molecules.  Again,  as  they  did  in  their 
uracil  study,  a refinement  of  42  force  constants  was  made  to  fit  ob- 
served frequencies,  and  the  average  frequency  error  ranged  from  12  - 
18  cm--*-.  They  pointed  out  that  an  empirical  assignment  for  methylated 
pyrimidine  derivatives  is  next  to  impossible  because  of  pronounced 
coupling  between  ring  stretching  modes  and  methyl  group  stretching  and 
deformation  modes. 

Detailed  studies  on  pyrimidine  bases  were  conducted  by  Tsuboi, 
Takahashi,  and  Harada  (1).  For  uracil,  cytosine,  adenine  and  guanine 
bases,  as  well  as  for  some  of  their  methylated  and  methylated-deute- 
rated  derivatives,  they  obtained  Raman,  preresonance  Raman  and  infrared 
spectra.  These  experimental  spectra  were  then  compared  with  those  of 
different  types  of  uridine  phosphoric  acids,  for  example,  8-uridine- 5 ' - 
phosphoric  acid.  Also,  they  carried  out  NCAs  for  these  molecules, 
including  JN  and  0 isotopic  shifts,  at  two  different  levels  of 
approximation  assuming  that  some  groups  behave  as  single  dynamical 
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units.  The  purpose  of  these  approximations  was  to  simplify  the  vibra- 
tional problem  to  one  with  less  atoms.  They  were  able  to  characterize 
Raman  lines  or  infrared  bands  describing  the  normal  modes  in  terms  of 
potential  energy  distributions.  The  force  constants  which  they  used  in 
their  calculation  were  obtained  by  a refinement  procedure. 

In  their  study  of  the  normal  modes  of  vibration  in  the  uracil 
residue  of  3-uridine-5 '-phosphoric  acid  Nishimura,  Hirakawa,  and 
Tsuboi  (6)  used  frequency  shifts  measured  for  a set  of  isotope 
derivatives  to  interpret  the  observed  vibrational  spectrum  of  uracil. 
They  found  that  the  N-H  in-plane  deformations  are  extremely  sensitive 
to  the  substitution.  The  agreement  between  the  observed  and 
calculated  15]sj_isotope  shifts  was  good  except  for  the  1397,  1689, 

1464,  and  1402  cm-'*'  resonance  Raman  lines  of  uracil.  They  explained 
that  this  disagreement  in  l^N-shifts  for  these  four  lines  could  be 
attributed  to  anharmonicity , Fermi  resonance,  and  vibrational  couplings 
with  the  solvent  molecules,  factors  which  were  not  taken  into  consider- 
ation in  their  calculation.  For  their  NCA  this  time  they  began  with 
force  constants  transferred  to  uracil  from  benzene. 

Very  recently,  a study  of  the  in-plane  vibrational  modes  in  uracil 
using  an  ab  initio  quantum  mechanical  calculation  was  made  by  Nishimura 
Tsuboi,  Kato,  and  Morokuma  (14).  A force  field  was  obtained  with  calcu 
lation  using  the  ST0-3G  basis  set,  which  was  then  scaled  to  that 
expected  for  a 4-31G  basis  set,  using  various  scaling  factors  estimated 
from  results  for  other  molecules  (acrolein,  urea,  and  formamide) . 

In  their  NCA  calculation  this  "corrected  to  4-31G  force  field"  could 
reproduce  observed  frequencies  of  uracil,  1 ,3-dideutero-,  5,6-dideu- 
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tero-,  and  perdeutero-uracil,  except  that  the  calculated  normal  fre- 
quencies are  8 % (on  the  average)  greater  than  the  observed  frequencies. 
However,  they  found  that  all  of  the  21  in-plane  calculated  frequencies 
were  readily  correlated  with  the  observed  frequencies,  and  that  the 
calculated  normal  modes  of  vibrations  were  not  contradictory  with  the 
available  experimental  data. 

In  order  to  study  the  band-splittings  in  the  double  bond  region  of 
the  spectra  of  uracil  and  methylated  uracils,  Bandekar  and  Zundel  (3) 
obtained  the  infrared  and  Raman  spectra  of  these  molecules  in  solution 
and  solid  phase.  Then  they  carried  out  a NCA  calculation  to  assign 
these  spectra,  using  a force  field  obtained  by  a least-squares  refine- 
ment procedure.  They  proposed  that  the  concept  of  transition  dipole- 
dipole  coupling  (3)  could  explain  the  observed  splitting  of  the  C=0 
double  bond  stretching  mode  in  these  molecules. 

Another  NCA  was  made  by  Bowman  and  Spiro  (15)  on  the  in-plane 
vibrations  of  uracil  and  three  isotopic  molecules  (N,N-dideutero-, 
C,C-dideutero-,  and  perdeutero-uracil)  using  a semi-empirical  molecular 
orbital  calculation  (MNDO)  to  obtain  a force  field  for  uracil.  This 
was  refined  to  reproduce  81  observed  frequencies  with  a mean  error  of 
9.7  cm  Then  they  applied  this  force  field  to  interpret  the  vibra- 
tional spectra  of  the  isoalloxazine  ring  of  flavin. 

Recently,  experimental  infrared  spectral  studies  of  uracil  and 
some  of  its  deuterated  derivatives  have  been  made  for  the  matrix- 
isolated  molecules  by  Maltese,  Passerini,  Nunziante-Cesaro,  Dobos  and 
Harsanyi  (2),  by  Barnes,  Stuckey,  and  Le  Gall  (5),  and  by  Szczesniak, 
Nowak,  Rostkowska,  Szczepaniak,  Person,  and  Shugar  (16).  Matrix 
isolation  spectroscopy  provides  a unique  method  for  study  of  vibra- 
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tional  spectra  of  monomeric  uracil  molecules  under  conditions  where 
interaction  between  molecules,  and  their  rotational  motions,  are  mini- 
mized (16-20).  This  technique  consists  in  trapping  the  uracil  molecule 
in  an  inert  solid  matrix  (argon,  xenon,  krypton  or  nitrogen)  at  10  K, 
and  obtaining  its  spectrum.  Raising  the  concentration  allows  the 
changes  in  the  spectrum  resulting  from  intermolecular  hydrogen  bonding 
to  be  observed  (20).  The  assignments  to  the  spectra  by  Maltese  et  al. 
(2),  Barnes  et  al.  (5)  and  Szczesniak  et  al.  (16)  will  be  considered 
in  detail  and  compared  in  Chapter  4 . 

Finally,  in  a recent  study  Lewis,  Miles,  and  Becker  (4)  obtained 

the  infrared  and  Raman  spectra  in  solid  phase  of  1-methyluracil  and 
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seven  isotopic  derivatives  ( H and  0) . In  an  attempt  to  distinguish 
in-plane  and  out-of-plane  modes  they  studied  polarized  Raman  and 
infrared  spectra  of  single  crystals.  Their  vibrational  assignments 
are  in  quite  good  agreement  with  previous  studies. 

Features  of  the  Present  Study 

In  this  study  we  shall  attempt  to  interpret  the  experimental 
matrix-isolated  spectra  of  uracils,  including  some  methylated  deriva- 
tives, by  carrying  out  NCA  calculations  using  force  constants  obtained 
at  the  ab  initio  level.  In  the  past,  as  described  in  the  previous 
section,  most  of  the  spectroscopic  studies  of  uracils  were  made  on 
solid  samples  or  solutions,  resulting  in  broad  infrared  spectra  in 
which  overlapping  bands  are  not  resolved,  due  to  intermolecular 
hydrogen  bonding,  interactions  with  polar  solvents,  etc.  The  advantage 
of  using  matrix-isolation  infrared  spectroscopic  data  is  that  the 
spectra  of  uracils  are  expected  in  general  to  be  equivalent  to  the  pure 
vibrational  spectra  for  the  isolated  monomeric  molecules;  hence  they 
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can  be  compared  directly  with  vibrational  spectra  predicted  theoreti- 
cally for  the  isolated  molecule.  Also,  the  use  of  complete  force 
fields  obtained  at  the  ab  initio  level  is  more  reliable  for  this  cal- 
culation of  theoretical  spectra.  This  procedure  is  different  from 
many  of  the  previous  studies  of  uracils,  where  valence  force  constants 
were  in  some  instances  arbitrary  values  forced  to  fit  the  experimental 
frequencies.  Evenmore,  very  few  interaction  force  constants  could  be 
evaluated  in  those  studies.  All  these  aspects  will  be  considered  in 
great  detail  in  Chapters  2 and  4. 

A very  important  aspect  of  this  dissertation  involves  the  transfer 
of  vibrational  parameters  — force  constants  and  atomic  polar  tensors 
(APTs) . The  APTs  are  intensity  parameters  which  relate  the  dipole 
moment  derivatives,  obtained  in  an  ab  initio  calculation,  and  the  abso- 
lute infrared  intensities,  as  will  be  discussed  in  Chapters  2 and  3. 

We  transfer  ab  initio  force  constants  and  APTs  obtained  for  uracil  and 
other  molecules  (acetone,  N-methylformamide,  and  N-methylacetamide)  to 
make  predictions  of  frequencies  and  intensities  of  the  methylated 
uracil  molecules.  This  calculation  will  be  carried  at  two  levels  of 
approximation:  one  in  which  the  methyl  group  is  treated  as  a single 
dynamical  unit,  and  other  in  which  it  is  properly  treated  as  a four- 
atom  unit.  We  will  compare  the  predicted  results  with  the  experimental 
infrared  spectra.  In  Chapter  5 we  will  examine  the  transfer  of  vibra- 
tional parameters,  and  in  Chapter  6 we  will  apply  the  transfer  of  these 
parameters  to  predict  the  spectra  of  the  methylated  uracils. 

Another  feature  of  the  present  study  is  that  we  predict  the  abso- 
lute infrared  intensities  in  the  infrared  spectra.  Then,  we  are  able 
to  compare  the  observed  spectra  of  the  uracils  with  both  calculated 
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frequencies  and  intensities.  Therefore,  we  can  be  confident  that  most 
of  our  assignments  are  correct.  An  illustration  of  the  calculations 
for  uracil  will  be  given  in  Chapters  3 and  4. 


CHAPTER  2 
TECHNIQUES 

This  chapter  deals  with  the  experimental  as  well  as  with  the  theo- 
retical methods  followed  in  the  present  study.  We  shall  first  describe 
the  experimental  procedures  that  were  used  to  obtain  the  infrared-active 
frequencies  and  intensities  of  matrix  isolated  uracil  and  its  deriva- 
tives. This  is  followed  by  a brief  description  of  the  calculation  for 
predicting  the  vibrational  spectra  of  the  same  molecules.  Since  these 
experimental  and  theoretical  techniques  have  been  discussed  in  detail 
elsewhere  (explicit  references  are  given  later  in  this  chapter)  our 
presentation  will  emphasize  only  those  considerations  that  apply  to 
the  present  work. 

Experimental  Studies  of  Infrared  Spectra 
of  Matrix-Isolated  Species 

The  experimental  matrix  isolation  infrared  spectroscopic  data  used 
in  this  study  were  obtained  by  K.  Szczepaniak  and  coworkers  (M.  J.  Nowak, 
M.  Szczesniak  and  H.  Rostkowska)  in  their  laboratory  at  the  Institute  of 
Physics  of  the  Polish  Academy  of  Sciences  in  Warsaw,  Poland.  The  experi- 
mental procedure  which  they  employed  to  obtain  the  observed  frequencies 
and  relative  intensities  of  the  infrared  spectrum  of  uracil  and  its 
methylated  derivatives  isolated  in  noble  gas  matrices  has  been  described 
in  Refs.  16  to  20. 

A brief  qualitative  summary  is  that  uracil  was  isolated  in  different 
matrices  (argon,  krypton,  xenon,  and  nitrogen) , and  methylated  uracils  in 
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argon  and  nitrogen,  at  matrix-to-solute  ratios  of  about  1000  : 1,  at  a 
temperature  of  approximately  10  K.  Then  the  infrared  spectra  of  these 
matrices  were  recorded  on  a Perkin-Elmer  580B  spectrometer  at  a reso- 
lution of  1 cm  1 . No  evidence  of  H-bonding  (for  example,  to  form 
dimers)  was  found  in  these  spectra,  indicating  that  the  species  present 
in  these  matrices  are  the  isolated  monomers. 

Calculation  of  Vibrational  Spectra 

In  order  to  calculate  the  vibrational  frequencies  of  a molecule, 

the  vibrational  eigenvalue  secular  equation  must  be  solved: 

GFL  = LA  . (2-1) 

Here  G is  the  inverse  kinetic  energy  matrix,  F the  force  constants 

matrix,  L the  tranformation  matrix  between  the  internal  coordinates  and 

the  normal  coordinates,  and  A the  eigenvalue  matrix  of  the  frequency 

2 2 2 

parameters  (\±  = 4ir  c v^,  where  c is  the  speed  of  light  constant  in 
cm  sec  1 and  the  wavenumber  in  cm  it  follows  that  v^(cm  ^)  = 

1303.1 /XT).  Here  the  secular  equation  (Eq.  (2-1))  has  been  set-up 
following  the  well-known  GF  matrix  method  by  Wilson,  Decius  and  Cross 
(21),  and  it  can  be  solved  by  methods  they  described.  Alternatively, 
Miyazawa  (22)  has  described  a simpler  method  to  solve  this  equation 
based  on  symmetry  considerations.  In  our  discussion  we  shall  follow 
Miyazawa 's  method. 

While  we  calculate  the  vibrational  frequencies  of  the  molecule 
under  study,  by  finding  the  eigenvalues  of  Eq.  (2-1),  it  is  desirable, 
as  well,  to  obtain  some  prediction  of  the  absolute  infrared  intensities 
for  each  fundamental  band.  Comparison  of  both  the  calculated  frequencies 
and  intensities  with  the  observed  spectrum  will  aid  in  the  correct 


assignment  of  the  latter. 
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The  atomic  polar  tensors  (APTs)  for  the  atoms  in  the  molecule  being 
studied  are  the  parameters  that  determine  the  infrared  intensities. 

These  APTs  are  the  derivatives  of  the  dipole  moment  with  respect  to 
displacements  of  the  atoms  in  a fixed-Cartesian  coordinate  system  (23) . 
In  this  section  we  shall  describe  the  use  of  the  APTs  to  calculate  the 
absolute  infrared  intensities  for  the  fundamental  vibrational  modes  of 
a molecule. 

Vibrational  Frequencies 

The  fundamental  frequencies  of  a molecule  are  calculated  by  a 
normal  coordinate  analysis  (NCA)  to  obtain  the  solution  of  the  secular 
Eq.  (2-1).  The  NCA  is  a rather  involved  subject  which  has  been  dis- 
cussed elsewhere  in  great  detail  (21,24-26),  so  here  only  the  general 
procedure  will  be  described.  For  a nonlinear  polyatomic  molecule  with 
N atoms  there  are  3N  - 6 normal  modes  of  vibration  (3N  - 5 if  the 
molecule  is  linear) . By  using  the  molecular  symmetry  it  is  possible  to 
identify  and  classify  these  modes  of  vibration  into  their  symmetry  group 
species  (24),  as  will  be  considered  below. 

It  will  be  convenient  to  point  out  several  important  coordinate 
transformations  that  appear  in  the  NCA.  These  have  been  summarized  in 
Table  2-1.  For  example,  it  is  shown  there  that  the  transformation  from 
space-fixed  Cartesian  displacement  coordinates,  X,  to  molecule-fixed 
internal  displacement  coordinates,  R,  is  given  by 

R = BX  (2-2) 

where  B is  a (3N  - 6)  x 3N  matrix.  Taking  advantage  of  any  symmetry 
that  the  molecule  may  have,  and  using  the  sets  of  R internal  coordinates 
as  basis  functions,  we  apply  the  projection  operator  technique  (27)  to 
set  up  linear  combinations  of  the  internal  coordinates  that  are  proper 


TABLE  2-1 

COORDINATE  TRANSFORMATIONS 
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symmetry  adapted  functions  transforming  according  to  the  appropriate 
irreducible  representations  of  the  symmetry  point  group  of  the  molecule. 
The  transformation  from  the  internal  coordinates  R to  the  symmetry 
coordinates  S is  then  expressed  as 

S = UR  (2-3) 

where  U is  a (3N  - 6)  x (3N  - 6)  orthogonal  matrix,  when  existing 
redundancies  have  been  eliminated.  The  redundancy  problem  will  be 
explicitly  dealt  with  when  we  describe  the  in-plane  symmetry  coordinates 
of  uracil  and  of  the  methylated  uracils.  In  terms  of  the  S symmetry 
coordinates,  the  Gs  and  Fs  matrices  of  the  vibrational  equation 
(Eq.  (2-1))  are  block  factored,  where  each  symmetry  block  is  associated 
only  with  the  subset  of  symmetry  coordinates  which  transform  according 
to  the  same  irreducible  representation  of  the  molecular  point  group. 

Another  important  relation  shown  in  Table  2-1  is 

Q = Ls_1S  (2-4) 

where  the  symmetry  coordinates  S are  transformed  to  normal  coordinates 
Q by  Ls  ^ , which  is  the  inverse  of  the  normal  coordinate  transformation 
matrix  Lg.  From  Eq.  (2-1)  it  can  be  seen  that  this  Ls  matrix  is  the 
eigenvector  matrix  obtained  from  the  solution  of  the  vibrational  secular 
equation. 

Overall  transformations  between  coordinates  are  also  given  in  Table 
2-1.  For  example,  the  transformation  from  Cartesian  coordinates  X to 
normal  coordinates  Q is 

Q = Ls_1UBX  = (2-5) 

where  the  matrices  in  the  first  part  have  been  defined  previously  in 
this  section,  and  Z’'  is  the  transpose  (equal  to  the  inverse  since  Z is 
orthogonal)  of  the  normal  coordinate  transformation  matrix  Z from 
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mass-weighted  Cartesian  coordinates  q to  normal  coordinates  Q,  and  If5  is 
a diagonal  matrix  whose  elements  are  the  square  roots  of  the  mass  of 
each  atom  (28) . Also  found  in  this  table  are  expressions  for  the 
kinetic  and  potential  energy  functions  in  terms  of  the  various  coordi- 
nate systems.  For  example,  the  relation 

2T  = X+MX  (2-6) 

gives  the  kinetic  energy  in  terms  of  the  time  derivatives  (velocities) 

of  the  Cartesian  coordinates  X,  the  diagonal  M matrix  defined  above, 

• • + 

and  the  transpose  of  X,  written  as  X . 

Two  important  matrices  in  the  NCA  are  the  G and  the  A matrices. 

The  first  one  is  defined  as 

G = BM_1B+  . (2-7) 

Here  B is  the  transformation  matrix  given  by  Eq.  (2-2),  B^  its  transpose 

and  M ^ the  inverse  of  the  diagonal  M matrix.  The  dimensions  for  the 
unsymmetrized  G matrix  are  (3N  - 6)  x (3N  - 6).  The  A matrix,  of 
dimensions  3N  x (3N  - 6),  transforms  the  internal  coordinates  R into  the 
Cartesian  coordinates  X 

X = AR  . (2-8) 

Substitution  of  R = BX  (Eq.  (2-2))  gives 

X = ABX  . (2-9) 

Therefore, 


AB  = I (2-10) 

where  I is  the  3N  x 3N  unit  matrix.  Premultiplying  G in  Eq.  (2-7)  by 
A gives 

AG  = ABM-1B+  = M-1B+  (2-11) 

or 


A = M 1B+G  1 


(2-12) 
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The  vibrational  secular  equation  in  terms  of  internal  coordinates 
(Eq.  (2-1))  is  most  conveniently  set-up  in  symmetric  form  according  to 
Miyazawa's  method  (22).  The  geometrical  basis  for  this  method  was 
discussed  by  Person  and  Crawford  (29).  According  to  this  method,  first 
consider  the  equation 

GY  = YT  (2-13) 

where  Y is  the  matrix  of  the  eigenvectors,  and  T is  the  diagonal 
eigenvalue  matrix,  of  G.  Since  G is  a real  symmetric  matrix,  Y is 
orthogonal;  thus 

G = YTY+  = = WW+  . (2-14) 

Here  we  have  defined  the  (3N  - 6)  x (3N  - 6)  matrix  W = YT^.  Using 
Eq.  (2-14)  the  secular  equation  in  Eq.  (2-1)  can  be  written  in  the 
form 

WW+FL  = LA  (2-15) 


or 

W+FW(W_1L)  = (W_1L)4 


If  we  define  a new  matrix  C as 

C = W_1L 


(2-16) 


(2-17) 


Eq.  (2-16)  takes  the  form 

W+FWC  = CA  . (2-18) 

Defining  also  H as 

H = W+FW  , (2-19) 


Eq.  (2-18)  becomes 

HC  = CA  . (2-20) 

Hence,  the  solution  of  the  secular  equation  (Eq.  (2-1))  is  simplified 
to  finding  an  orthogonal  matrix  C which  diagonalizes  the  symmetric  H 
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matrix,  since  it  follows  from  Eq.  (2-20)  that 


L 1 = A -1L  +F 
~ s ~ s - s ~ s 


A = CTHC  . (2-21) 

The  solution  of  Eq.  (2-1)  is  thus  reduced  to  diagonalization  first  of 

the  symmetric  G matrix  (Eq.  (2-14)),  calculating  W and  then  H,  and 

diagonalizing  that  symmetric  matrix.  Diagonalization  methods  for 

symmetrical  matrices,  based  on  iterative  methods  like  Jacobi's  (30), 

are  very  well  suited  for  computer  programs. 

It  is  interesting  to  note  that  once  the  secular  equation  in 

symmetric  form  has  been  solved,  the  Ls  transformation  matrix  can  be 

obtained  from  Eq.  (2-17) 

L = W C = Y T . (2-22) 

~s  ~s~s  ~s~s  ~ s 

And  from  Eqs.  (2-18)  and  (2-22) 

(2-23) 

The  units  that  we  use  in  the  NCA  deserve  some  discussion.  The 
Cartesian  displacement  coordinates  of  the  atoms  are  expressed  in 

angstroms  (&,  1 & = 10  ® cm) , so  the  internal  displacement  coordinates 

o 

for  bond  stretches  (Ar)  have  units  of  A and  the  internal  coordinates 

for  angle  bends  (A9)  are  in  radians.  Hence,  the  B matrix  element  for 

a bond  stretch  is  dimensionless,  while  that  for  a bending  internal 

0_1 

coordinate  will  have  units  of  A . Since  we  use  atomic  mass  units 

• —27 

(u,  1 u = 1.66056  x 10  kg)  for  the  masses,  the  G matrix  elements  all 

-1  o—i  o_2 

have  units  of  u . The  units  of  the  F matrix  are  mdyne  A (=  aJ  A ) 

o 

for  the  force  constants  involving  two  stretches,  mdyne  A for  bend-bend 
elements,  and  mdyne  for  bend-stretch  elements. 

Potential  Energy  Distributions 

After  having  carried  out  the  NCA,  it  is  useful  to  calculate  what 
percentage  of  the  potential  energy  associated  with  a given  normal 
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coordinate  is  contributed  by  each  symmetry  coordinate.  From  Table  2-1, 
the  potential  energy  2V  is  given  by 

2V  = Q+AQ  . (2-24) 

Since  the  eigenvalue  matrix  A is  given  by 

A = L+FL  , (2-25) 

substitution  of  this  relation  into  Eq.  (2-24)  gives 

2V  = Q+LfFLQ  . (2-26) 

For  the  kth  normal  coordinate,  whose  frequency  parameter  is  Xk, 

\ - 2Hi+FijLjk " 2Fi3LikLjk  • <2-27) 

i.j  i.j 

The  total  potential  energy  distribution  (PED)  for  the  kth  normal  coordi- 
nate was  defined  by  Morino  and  Kuchitsu  (31)  as 

V(Qk)  = (^)Qk2ZFijLikLjk  . (2-28) 

i.j 

using  Eqs . (2-24)  and  (2-27).  Examining  the  magnitudes  of  the  terms 

contributing  to  the  sum  in  Eq.  (2-28),  it  is  argued  that  in  general, 

2 

the  value  of  F^jL^Lj^  is  large  only  when  i = j.  Thus,  the  F^L-j^ 
terms  are  usually  the  most  important  in  determining  the  PEDs . There- 
fore, it  has  become  conventional  to  define  the  percentage  PED  from  the 

2 

ith  symmetry  coordinate  to  be  given  by  the  ratio  of  F^L^k  to  the  total 

sum  T F Ljj.Lj k (this  sum  is  equal  to  Ak)  times  100 

i.j  o 

PED±k  = (Fi;LLik  /\)  x 100  . (2-29) 

Here,  F^  is  the  diagonal  force  field  element  associated  with  the  ith 

symmetry  coordinate,  is  the  eigenvalue  of  the  kth  normal  coordinate, 

and  Lj[k  is  an  element  from  the  normal  coordinate  transformation  matrix  L. 

The  PEDs  provide  information  about  the  contribution  of  each 

symmetry  coordinate  ± to  the  potential  energy  for  the  kth  normal  mode  of 

vibration.  With  the  aid  of  Eq.  (2-29) , we  can  determine  whether  all  the 

symmetry  coordinates  contribute  equally  to  the  total  potential  energy  of 
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a given  normal  coordinate  or  whether  just  a few  predominate.  It  is  a 
valuable  guide  to  the  interpretation  of  results  from  the  NCA. 

Infrared  Intensities 

The  other  aspect  of  the  vibrational  spectra  that  we  want  to 
consider  is  the  absolute  infrared  intensity  of  each  absorption  band. 

For  a normal  vibrational  mode  k the  absolute  infrared  intensity  is  pro- 
portional to  the  square  of  the  molecular  dipole  moment  derivative  with 
respect  to  the  kth  normal  coordinate  (23)  . The  components  along  the 
space-fixed  Cartesian  axes  of  the  dipole  moment  derivatives  with  respect 
to  the  normal  coordinates  can  be  written  as  a 3 x (3N  - 6)  tensor, 
designated  as  Pq.  In  order  to  interpret  or  to  predict  these  values,  the 
first  thing  to  do  is  to  calculate  the  APTs  Px  in  the  space-fixed  Carte- 
sian coordinate  system.  In  the  next  section  we  will  describe  the  method 
used  to  calculate  Px.  Let  us  now  describe  the  relationship  between 
these  tensors  and  the  absolute  infrared  intensities. 

The  relation  between  Pq  and  Px  is  given  by  (32) 

PQ  = PXAL  (2-30) 

where  the  A and  L matrices  are  obtained  through  the  NCA,  as  previously 
described.  The  elements  in  Px  have  units  of  electrons  (e) , those  in  A 
are  unitless  and  those  in  L have  units  of  u . The  resulting  units  of 

-h 

Pq  are  eu 

For  absorption  of  light  by  a molecule  during  a transition  from  the 
± to  the  i + 1 vibrational  state  of  the  sth  normal  mode,  the  absolute 
absorption  coefficient  Ag  is  given  by  (33) 

gSj 

As  = (8Tr2N0vs/3hc)  2 1 + 1 1 P I i>  I 2 (2-31) 

where  NQ  is  Avogadro's  number,  vs  is  the  transition  frequency  in  cm 

-1 

h is  Planck's  constant,  c is  the  speed  of  light  in  cm  sec  , p is  the 
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dipole  moment  operator,  and  <i  + l|p|i>  is  the  transition  moment 
integral.  The  sum  on  1 is  over  the  gg  degenerate  states;  all  of  the 
vibrational  modes  in  the  uracils  are  singly  degenerate  (gs  = 1) . 

Under  the  assumption  that  mechanical  anharmonicity  is  small,  the 
wavefunctions  are  perfect  harmonic  oscillator  functions.  After 
expanding  the  dipole  moment  in  a Taylor ' s series  in  the  normal 
coordinates  Q 

p = p°  + S(3p/9QJQ-  + ^ 2T  Il02p/9Q  3Q.)Q  Q.  + ...  (2-32) 

s ss  sj  b J b J 

the  second  and  higher  order  terms  are  neglected  (assumption  that 
there  is  no  electrical  anharmonicity) , so 

As  = (N0Trgs/3c2)|3p/9Qs|2  • (2-33) 

Here  |9p/9Qg|2  is  given  by 

| 9p/9Qs | 2 = (9px/3Qs)2  + (3py/9Qs)2  + (9pz/9Qg)2  (2-34) 

where  |9p/9Qg|  is  the  absolute  value  of  the  dipole  moment  vector 
derivative  with  respect  to  the  sth  normal  coordinate,  and  3px/9Qg, 
9py/9Qg  and  9pz/9Qg  are  the  components  of  9p/9Qg  along  each  of  the 
x,  y and  z space-fixed  Cartesian  axes. 

Evaluation  of  the  constants  in  Eq.  (2-33)  when  Ag  is  in  km  mol  ^ 
and  9p/3Qg  in  eu-15  units  gives  (32) 

Ag(km  mol"1)  = 974.9  gs|9p/9Qs|2  . (2-35) 

Now  the  calculation  of  the  absolute  infrared  intensities  is 
straightforward.  The  calculated  APTs  in  space-fixed  Cartesian  system 
(Px)  are  transformed  into  normal  coordinate  space  (Pq)  with  the  use  of 
Eq.  (2-30),  where  the  A and  L matrices  have  been  obtained  from  the 
NCA.  Then  the  Pq  values  are  substituted  in  Eq.  (2-35)  to  get  the  A 
values  in  km  mol  ^ . 
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Calculation  of  Atomic  Polar  Tensors 
The  atomic  polar  tensors  (APTs)  for  uracil  were  calculated  by 
S.  Chin  using  ab  init io  quantum  mechanical  methods  (17,18).  To 
summarize  briefly,  what  was  done  was  to  calculate  the  wavefunctions , 
energy,  and  dipole  moment  components  of  uracil  at  the  experimental 
equilibrium  configuration  (17)  and  then  to  repeat  the  calculation  at  the 
distorted  geometries  to  obtain  the  dipole  moment  derivatives  P^  by  the 
method  of  finite  differences  (34) . These  dipole  derivatives  are  calcu- 
lated with  respect  to  a set  of  mass-weighted,  molecule-fixed,  principal 
Cartesian  displacement  coordinates  5 orthogonal  to  the  translations  and 
rotations  and  reflecting  the  molecular  symmetry  (28).  The  orthogonality 
of  the  S coordinates  is  ensured  by  defining  them  to  be  orthogonal  to 
each  other  and  to  the  six  Eckart  conditions  (35,36),  which  define  three 
translations  Pt  of  the  entire  molecule,  and  three  rotations  Pr. 

The  APTs  in  the  principal  Cartesian  coordinate  system  are  given 
by  the  following  equation  (28,37) 

P = (P  iPcOUir2  (2-36) 

where  P is  a 3 x 6 matrix  that  consists  of  the  translational  and 
~P 

rotational  corrections  to  the  dipole  moment  derivatives  (38) , and  is 
equal  to 


1° 

0 

0 

0 

n°  /I  ^ 
P z'  J-yy 

-py^' 

?p  = 

0 

0 

0 

-PVlxx15 

0 

PVIzz1S 

(2-37) 

1° 

0 

0 

pyw5 

-pVV5 

0 } 

• 

Here,  the  p°s  are  the  components  of  the  permanent  dipole  moment  along 
the  principal  Cartesian  coordinate  axes  at  equilibrium  configuration; 
and  Ixx,  Iyy  and  Izz  are  the  principal  moments  of  inertia.  Also,  in 
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Eq.  (2-36)  P refers  to  the  dipole  derivatives  with  respect  to  the  S 
coordinates,  defined  as  (28) 


/; 


3pV  353n\ 


(2-38) 


3PX/3S7  3Px/3Sg  3px/3S9  .. 

9py/9S7  9Py/9Sg  9py / 95g  ...  9py/9S3N 

W3S7  3p2/9S8  8pz/9S9  ...  3p  z/3S3[|J 

The  U matrix  in  Eq.  (2-36)  is  the  3N  x 3N  orthogonal  transformation 
matrix  that  transforms  the  mass-weighted  principal  Cartesian  coordi- 
nates g into  the  5 and  P coordinates,  according  to 


6) 


s = u ? 


The  M matrix  has  been  defined  previously. 


(2-39) 


Finally,  the  full  polar  tensor  P calculated  by  Eq.  (2-36)  is 
given  by  this  general  form  (23,28,39) 


K/9X1 

V3T1 

9px/3zi 

9px/9X2 

3p  / 3Y9  ... 
rx  2 

9px/9Z 

p = 
~x 

9p  / 3X, 
ry  1 

9py/9Z1 

9Py/9X2 

3Py / 9Y2  ... 

9py/3Z 

lV3xi 

8pz/8Yj 

3pz/3zl 

9pz/9X2 

3Pz/3y2  • • • 

3pz/3Z 

(2-40) 


In  this  arrangement  the  first  three  columns  contain  the  APT  of  the 
first  atom,  the  next  three  the  APT  of  the  second,  and  so  on. 

A full  account  on  the  calculation  of  APTs  can  be  found  in 
references  38  and  39.  We  will  describe  some  properties  of  APTs  for 
uracil  in  Chapter  3. 

Transfer  of  Parameters  to  Predict  Infrared  Spectra 
Since  we  are  interested  in  making  quantitative  predictions  of  the 


infrared  spectra  of  the  polyatomic  methylated  uracil  molecules,  we  are 
particularly  interested  in  the  possibility  of  transferring  vibrational 
parameters  from  similar  smaller  molecules  or  molecular  fragments  to 
succesfully  achieve  this  purpose.  The  molecular  parameters  that  may  be 
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transferred  are  the  force  constants  and  the  APTs.  In  the  following 
discussion  we  shall  describe  how  this  technique  was  applied  to  the 
1-methyluracil  (1-Me  U)  molecule,  as  an  illustration.  First  we  will 
address  the  problem  of  the  transfer  of  force  constants,  and  later  the 
transfer  of  APTs. 

Transfer  of  ab  initio  force  constants 

Let  us  first  consider  a simple  scheme  to  show  how  the  transfer  of 
vibrational  parameters  should  work.  In  principle,  we  want  to  transfer 
the  parameters  from  some  molecules  with  similar  structural  fragments 
to  the  molecule  under  study.  For  the  specific  case  of  1-Me  U,  Fig.  2-1 
depicts  in  a schematic  way  the  strategy  to  be  followed  in  the  transfer 
of  parameters  when  the  Hy  atom  in  uracil  is  replaced  by  a methyl  group 
to  form  1-Me  U. 

As  seen  in  Fig.  2-1,  two  simple  approaches,  indicated  as  Calcula- 
tion I and  Calculation  II,  can  be  followed  in  the  procedure.  In 
Calculation  I the  methyl  group  is  considered  to  be  just  a single 
dynamical  unit  of  mass  15  u.  In  other  words,  1-Me  U is  regarded  as  a 
12  atom  molecule,  similar  to  the  calculation  of  1-deuterouracil  (17). 

In  Calculation  II  the  methyl  group  is  properly  treated  as  a four-atom 
CH^  unit,  and  the  predictions  for  1-Me  U include  the  vibrations  of  the 
methyl  group  as  well  as  those  of  the  ring. 

The  force  constants  that  we  used  in  Calculation  I,  except  for  the 
N^-Me  values,  were  transferred  from  the  ab  initio  calculation  with  the 
ST0-3G  basis  set  "converted"  to  values  expected  at  the  4-31G  level  for 
uracil  obtained  (and  described)  by  Nishimura,  Tsuboi,  Kato,  and 
Morokuma  (14).  Of  course,  we  expect  that  the  force  constants  for  the 
bonds  in  both  rings  (in  uracil  and  1-Me  U)  will  have  nearly  the  same 
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Figure  2-1.  Scheme  for  Calculations  I and  II  of  1-Methyluracil  (1-Me  U) . 

Vibrational  parameters  (ab  initio  force  constants  and  APTs) 
are  transferred  to  1-Me  U from  uracil,  N-Methylformamide, 
N-Methylacetamide  and/or  acetone. 
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values,  given  their  structural  similarity.  We  assumed  also  that  the 
structural  parameters  (bond  distances  and  angles)  of  both  rings  are  the 
same.  For  the  Nj-Me  bond  the  values  for  the  stretching  and  bending 
force  constants  were  transferred  directly  from  the  N^-Cy  bond  of 
N-methylformamide,  from  the  ab  initio  calculation  by  Sugawara,  Hirakawa, 
Tsuboi,  Kato,  and  Morokuma  (40).  Also,  we  used  the  Nj-Me  bond  length 
that  they  reported  as  the  experimental  value.  The  interaction  force 
constant  values  for  the  Nj-Me  stretch  with  the  ring  stretches  in  1-Me  U 
were  assumed  to  be  the  same  as  those  between  Nj-Cy  and  the  OCN  group 
in  N-methylformamide.  In  Chapter  4 we  present  the  "ab  initio  ST0-3G 
converted  to  4-31G"  force  constants  for  uracil  obtained  by  Nishimura 
et  al.  (14),  and  in  Chapter  6 all  the  force  constants  (transferred  from 
uracil  and  N-methylformamide)  for  the  methylated  uracils  will  be  given. 

In  Calculation  II  the  force  constants  for  the  pyrimidine  ring  of 
1-Me  U were  again  taken  from  the  ab  initio  calculation  by  Nishimura 
et  al.  (14).  For  the  sake  of  consistency,  we  decided  to  use  force 
constants  for  the  N^-CH^  fragment  obtained  by  ab  initio  methods.  For 
this  purpose,  the  values  that  Sugawara  et  al.  (40)  obtained  for 
N-methylformamide  at  the  4-31G  level  seemed  very  appropriate  to  use. 
Therefore,  the  set  of  force  constants  for  the  N^-CH^  group  in  N-methyl- 
formamide  from  their  ab  initio  calculation  (40)  were  transferred  to 
the  Nj-CH^  group  in  1-Me  U. 

The  procedure  just  described  was  also  used  for  3-methyluracil 
and  1 , 3-dimethyluracil . All  the  resulting  force  fields  are  given 


in  Chapter  6. 
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Transfer  of  APTs 

For  Calculation  I we  transferred  directly  all  the  atomic  polar 
tensors  from  uracil  to  1-Me  U,  including  the  APT  for  the  Me  unit,  which 
is  that  for  H7  of  uracil,  so  that  the  APTs  satisfy  the  null  condition 
(the  sum  of  all  the  APTs  of  a molecule  is  the  null  tensor) . Therefore, 
the  problem  of  predicting  intensities  for  1-Me  U in  Calculation  I is 
similar  to  that  of  predicting  the  intensities  for  one  of  the  deuterated 
derivatives  of  uracil  (17). 

It  should  be  pointed  out  that  the  assumption  that  the  APT  for  H7 
atom  of  uracil  is  the  same  as  that  for  the  methyl  group  is  not  very 
realistic,  so  that  Calculation  I is  expected  to  give  only  approximate 
results,  for  predicted  intensities  as  well  as  for.  frequencies. 

For  Calculation  II  the  H7  atom  in  uracil  is  replaced  by  a CH3  group 
to  form  1-Me  U.  In  doing  this  we  considered  that  the  APTs  for  the 
methyl  ( — CH3)  fragment  could  be  taken  from  the  atoms  in  the  methyl 
group  attached  to  the  nitrogen  atom  in  N-methylacetamide.  Among  other 
things,  the  CH3  group  in  both  molecules  (1-Me  U and  N-methylacetamide) 
is  attached  to  a nitrogen  atom,  with  a C=0  group  as  a neighbor  bond. 
Also,  the  well  known  conjugate  character  of  the  C7-N3  bond  in  N-methyl- 
acetamide  could  be  expected  to  give  almost  the  same  rigidity  to  that 
molecule  as  is  found  in  the  N1-C7  bond  of  1-Me  U (see  scheme  in  Fig.  2-1 
for  the  numbering  of  atoms).  Accordingly,  we  decided  to  undertake  a 
separate  study  of  N-methylacetamide  in  order  to  obtain  its  APTs,  and 
then  transfer  the  APTs  for  the  N-CH3  group  of  this  molecule  to  the 
N-CH3  group  of  1-Me  U.  We  present  the  APTs  for  N-methylacetamide 
in  Chapter  5. 
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As  is  also  illustrated  in  the  scheme  in  Fig.  2-1,  we  also  tested 
the  sensitivity  of  the  predicted  intensities  to  the  APIs  by  calculating 
the  intensities  using  the  APTs  of  the  C-CH^  group  of  N-methylacetamide 
for  the  intensity  parameters  assumed  for  the  N^-CH^  group  of  1-Me  U. 
Another  calculation  was  done  using  the  APTs  transferred  from  the 
C-CH^  groups  of  acetone,  taken  from  Rogers  (39). 

Because  the  null  condition  must  be  satisfied  in  all  the  intensity 
predictions  of  Calculation  II  described  above,  we  chose  the  APT  for 
either  or  Cy  in  order  to  satisfy  it,  transferring  all  other  APTs 
from  uracil,  N-methylacetamide  or  acetone,  and  then  compared  the 
calculated  intensities  with  the  observed  spectrum.  Applying  the 
appropriate  rotation  of  the  coordinates  to  these  APTs  made  them 
transferable  to  3-Me  U and  1,3-di  Me  U.  The  different  sets  of  APTs 
that  were  obtained  for  all  these  molecules  will  be  given  in  Chapter  6. 


CHAPTER  3 

ATOMIC  POLAR  TENSORS 

This  chapter  deals  with  the  atomic  polar  tensors  (APTs)  for  the 
uracil  molecule  and  their  use  in  the  prediction  of  absolute  infrared 
intensities  for  this  molecule.  Before  considering  any  specific  data, 
we  shall  review  some  additional  aspects  of  atomic  polar  tensor  theory. 
Then  we  shall  present  the  calculated  APTs  for  uracil,  and  review  the 
predicted  absolute  infrared  intensities  for  uracil  and  examine  some 
sum  rules.  In  particular  we  shall  examine  separately  the  total 
intensity  for  the  in-plane  modes  and  that  for  the  out-of-plane 
vibrational  modes,  compared  with  the  calculated  ratio  and  with  the 
experimental  value  of  this  ratio. 

Further  Theory  of  Infrared  Intensities 
As  discussed  in  Chapter  2,  the  set  of  normal  coordinates  des- 
cribes the  vibrational  modes,  and  can  be  expressed  in  terms  of 
symmetry  coordinates  (38)  as 

9p/9Qa  = ^-,LL6  (3p/3Sj)  (3-1) 

4 i=l 

or  conversely 

9p/9Si  = (9p/9Q*)  (3-2) 

i=l 

where  LL6  is  the  i,^th  element  of  the  L matrix,  and  Lii-^  is  the  a ,ith 

element  of  the  inverse  matrix  of  L. 

Since  we  cannot  determine  the  sign  of  9p/9Q/1  by  experiment,  the 

use  of  Eq.  (3-2)  leads  to  more  than  one  possible  set  of  values  for  the 

-+  k-1 

dipole  derivatives  9p/9Si  (in  fact  to  2 values,  where  k is  the 
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number  of  terms  (3N  - 6)  in  the  summation).  The  summation  in  Eq.  (3-2) 
is  made  separately  for  all  modes  of  a given  symmetry;  for  a symmetry 
block  with  k modes,  there  will  be  2 possible  sets  of  8p/9S^  values. 

In  order  to  obtain  information  for  the  determination  of  the  correct 
signs  for  these  values,  intensity  data  from  isotopically  substituted 
molecules  are  usually  used  (41).  Also,  the  analysis  of  the  Coriolis 
interactions  of  the  molecule  sometimes  offers  valuable  information 
about  such  signs  (42,43).  Different  research  groups  have  demonstrated 
that  one  possible  method  of  overcoming  this  signs  problem  lies  in  the 
use  of  quantum-mechanical  calculations  (38,44-47).  In  this 
dissertation  we  are  using  calculated  APTs,  with  the  signs  from  their 
calculation,  to  predict  the  derivatives  in  Eq.  (3-1),  so  we  do  not 
have  this  sign  problem,  which  arises  when  we  use  experimental  intensi- 
ties to  derive  values  of  9p/9S^. 

A number  of  theories  have  been  proposed  for  expressing  3p/9Q^  in 
terms  of  parameters  characterizing  the  individual  bonds  or  atoms  of  the 
molecule.  A great  deal  of  interest  during  the  last  decade  has 
concerned  the  development  and  application  of  such  approaches  to  molecules 
in  the  gas  phase  (48).  These  include  the  valence-optical  theory, 
developed  by  Gribov  (49),  Sverdlov  et  al.  (50),  Gussoni  (51),  and 
other  workers  (52-54).  Here,  the  dipole  moment  is  expressed  as  a sum 
of  bond  dipole  moments  all  directed  along  their  respective  bonds.  The 
components  of  these  bond  dipole  moments  and  their  derivatives  with 
respect  to  the  internal  coordinates  are  called  electrooptical  parameters 
(eop's).  The  major  problem  with  this  approach  is  that  it  contains  too 
many  parameters  (greater  than  the  number  of  independent  equations)  to  be 
evaluated  from  experiment,  when  used  in  its  purest  form.  It  can  be 
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simplified,  however,  to  give  considerable  insight  for  the  interpretation 
of  intensities  (55).  For  a review  on  eop's  see  Ref.  51. 

Some  modifications  to  the  valence-optical  theory  have  been 
suggested  (56-58).  For  example,  in  the  hybrid  orbital  rehybridization 
model  of  Wiberg  and  Wendoloski  (56)  the  9p/SQA  values  are  described  in 
terms  of  a contribution  resulting  from  a geometric  reorientation  of  the 
bond  dipole  moment,  and  a second  term  which  describes  the  change  in 
bond  dipole  moment  due  to  orbital  rehybridization.  All  the  parameters 
involved  in  such  terms  are  expressed  as  simple  functions  of  the  bond 
hybridization. 

Other  similar  theories  based  on  effective  charges  of  bonds  or 
atoms  have  been  proposed  (57,58).  The  bond  dipole  moment,  as  defined 
in  the  valence-optical  theory,  can  be  expressed  as  a product  between 
an  effective  bond  charge  and  the  bond  vector.  Then,  differentiation  of 
the  bond  dipole  moment  with  respect  to  the  internal  coordinates  gives  a 
bond  charge  derivative.  This  theory  has  been  extensively  developed  and 
applied  by  Smit  and  van  Straten  (58) . 

Another  intensity  theory  which  has  had  considerable  impact  during 
recent  years  is  the  Atomic  Polar  Tensor  theory  (APT) . It  was  suggested 
by  Biarge,  Herranz  and  Morcillo  in  1961  (59),  and  reviewed,  reformulated 
and  extensively  applied  by  Person  and  Newton  (23,32,60).  Reviews  on  the 
formulism  and  applications  of  this  theory  are  found  in  Refs.  38,  48,  and 
54,  and  we  shall  summarize  some  aspects  in  the  following  discussion. 

In  the  APT  theory  the  components  of  the  dipole  moment  derivatives 
are  evaluated  with  respect  to  a space-fixed  Cartesian  coordinate  system 
(Eq.  (2-40)).  An  advantage  with  APTs  is  that,  because  they  are 
expressed  in  space-fixed  coordinates,  they  are  now  invariant  to 
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isotopic  substitution,  so  we  can  expect  the  APTs  to  transfer  from  one 
isotopic  form  of  a molecule  to  another  (38).  Also,  the  APTs  are 
calculated  directly  in  the  quantum  mechanical  calculations,  which  are 
usually  made  in  a space-fixed  coordinate  system.  A distinction  between 
APTs  and  eop’s  is  that  APTs  are  "intensity  parameters"  that  have  a 
one-to-one  relationship  to  the  experimental  intensity  values.  On  the 
other  hand,  the  eop's  are  useful  to  interpret  these  intensity  parame- 
ters. Using  group  theory,  Decius  and  Mast  (57)  proved  that  the  number 
of  independent  parameters  in  the  APTs  for  a molecule  is  the  same  as 
the  number  of  equations  used  to  determine  them  from  the  experimental 
data.  This  aspect  of  APTs  has  been  discussed  by  Person  (38)  in  more 
detail. 

The  APTs  for  each  atom  a (Eq.  (2-40)),  in  a molecule  composed  of 

N atoms,  are  collected  in  the  full  3 x 3N  polar  tensor  P 

~x 


The  expression  of  the  APTs  in  space-fixed  coordinates  (Px)  is  particu- 
larly useful  because  it  is  possible  to  transform  them  into  normal 
coordinate  space  (Pq)  according  to  Eq.  (2-30),  and  thus  to  obtain  the 
predicted  absolute  intensities  from  P^. 

The  APTs  are  always  expressed  relative  to  a particular  Cartesian 
coordinate  system.  If  the  APT  for  an  atom  in  one  molecule  is  to  be 
compared  or  transferred  to  the  same  kind  of  atom  in  another  molecule, 
it  is  necessary  to  express  both  APTs  in  the  same  coordinate  system. 

The  "bond  coordinate  system"  has  proven  to  be  a useful  system  for  such 
purposes  (60).  In  this  system  we  define  the  positive  x,  y,  or  z direc- 
tion to  be  along  the  bond  in  the  direction  from  the  central  to  the 


terminal  atom  for  the  two  bonds  in  consideration. 
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In  general,  by  applying  an  orthogonal  rotation  transformation  T 
the  "bond  coordinate  axes"  can  be  rotated  with  respect  to  the  "molecular 
coordinate  axes";  hence  the  APTs  in  one  axes  system  are  related  to  the 
APTs  in  the  other  system  by  (28) 

Pxa(new)  = T{Px“(old) }T+  . (3-4) 

Here  T is  a transformation  matrix  (T^  is  its  transpose)  whose  elements 
are  the  cosines  between  the  old  and  new  axes. 

In  addition,  the  APTs  for  the  N atoms  in  a molecule  are  also 

related  through  the  null  condition 

zk.  Px  = 0 (3-5) 

a=i 

which  arises  because  the  pure  translation  of  a molecule  does  not  alter 
the  dipole  moment  vector  of  the  molecule  and  hence  a vanishing  sum 
results  (23). 

The  invariant  APT  properties  (32,28,61)  associated  with  absolute 
infrared  intensities  are  the  effective  atomic  charge  xa>  the  mean 
dipole  moment  derivative  pa,  and  the  atomic  anisotropy  8a  . An  APT 
invariant  is  independent  of  the  choice  of  the  Cartesian  coordinate 
system  used.  The  effective  charge  Xa  is  related  to  the  APT  by  the 
following  relation  (32,61) 

Xa2  = (l/3)Tr{Pxa(PxV}  (3-6) 

where  Pxa  is  the  APT  for  atom  a,  (Pxa)^  its  transposed  form,  and  Tr(A} 
designates  the  trace  of  the  A tensor.  Actually,  the  right -hand-side 
of  Eq.  (3-6)  is  just  one- third  of  the  sum  of  the  squares  of  all 
elements  of  the  Pxa  tensor  for  the  a atom. 

The  mean  dipole  moment  derivative  pa  is  given  by  (28,61) 

Pa  = (l/3)Tr{Px“)  . (3-7) 

Also,  the  anisotropy  8a^,  which  measures  the  deviation  of  Px  from  a 


36 


constant  diagonal  tensor  of  value  pa  is  given  by  the  following 


relationship  (61) 


where,  for  example,  pxz  stands  for  3px/9Z  . The  relation  between  xa> 


It  has  been  found  that  Xa  i-s  nearly  the  same  for  all  hydrogen 
atoms  of  a number  of  hydrocarbons,  with  some  exceptions  (61).  Indeed, 
King,  Mast  and  Blanchette  (62)  suggested  that  the  effective  charge  for 
an  atom  may  be  a transferable  intensity  parameter.  Subsequent  work 
(32,45,55)  suggests  this  idea  may  be  somewhat  simplified,  but  it  may 
be  true  for  atoms  in  similar  chemical  environments.  In  a comparative 
study  of  effective  atomic  charges  for  hydrogen,  fluorine  and  oxygen 
atoms  in  a series  of  compounds,  King  (61)  concluded  that  Xa  values  for 
these  atoms  in  such  diverse  compounds  do  fall  within  a narrow  ranges 
of  values.  These  ranges  are  Xy  = 0.101  ± 0.005  e in  all  hydrocarbons 
except  HCCH,  XH  = 0.072  ± 0.012  e in  the  methyl  halides,  XF  = 0.595  ± 
0.039  e,  and  X0  = 0.624  ± 0.099  e.  He  pointed  out  that  these  results 
suggest  that  xa  represents  a particularly  localized  combination  of 
chemical  bond  properties. 

A useful  expression  relating  the  xa  values  for  the  atoms  in  a 
molecule  to  the  sum  of  the  absolute  intensities  for  all  the  vibrational 
modes  of  that  molecule,  was  derived  by  King,  Mast  and  Blanchette  (62) 


pa  and  6a2  is  (61) 


Xa2=  (Fa)2+  (2/9)  (Ba)  2 


(3-9) 


(3-10) 


Here,  m is  the  mass  of  atom  a,  and  £1  is  the  rotational  correction  (62) 
a 
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where  Pp  was  previously  defined  in  Eq.  (2-37).  Evaluating  Eq.  (3-11), 
we  find 

a - {<Pz°>2/Ixx  + <Pz°>2/Iyy  + l<Px°>2  + <Py°)2]/Iz2} 

+ {(py°)2/Ixx  + (px°)2/Iyy}  • (3-12) 

Each  one  of  these  quantities  was  defined  earlier  when  we  introduced 
Eq.  (2-37). 

In  Eq.  (3-10)  the  values  of  Xa^/ma  are  usually  in  e^u-1  units  when 
the  APT  values  are  expressed  in  e units,  and  masses  in  atomic  mass 
units  (u) . Hence,  the  terms  in  ft  must  be  evaluated  using  dipole 
moments  in  eA  (le&  = 4.8  Debye)  and  moments  of  inertia  in  u^ . When 
As  is  in  units  of  km  mol-^  King's  intensity  sum  rule  (Eq.  (3-10)) 
becomes  (32) 

2^TAS  (km  mol-*)  = 974 .9^51 3xa^/%  - £2^  • (3-13) 

We  can  use  a relationship  similar  to  Eq.  (3-13)  for  the  sum  of  the 
infrared  intensities  of  modes  belonging  to  each  species  contained  in 
the  symmetry  group  to  which  the  molecule  belongs  (63,64).  We  shall 
discuss  this  topic  in  more  detail  after  presenting  the  APTs  for  uracil. 

Finally,  we  should  consider  the  quantum-mechanical  description  of 
the  APTs  given  by  King  and  Mast  (65)  . They  interpreted  the  APTs  as  the 
sum  of  three  contributions:  a net  charge  (calculated  as  the  Mulliken 
charge),  the  charge  flux,  and  the  "interference",  which  is  now  called 
the  "overlap"  term  (66) . This  treatment  has  been  slightly  modified 
by  Zilles  (66)  to  obtain  this  expression 

p a _ p “(charge)  + P “(charge  flux)  + P “(overlap)  . (3-14) 

The  physical  significance  of  Eq.  (3-14)  is  the  following:  Px“(charge) 
is  the  contribution  to  the  dipole  derivative  due  to  motion  of  the 
fixed  equilibrium  charge,  Pxa (charge  flux)  is  due  to  the  change  in 
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charge  on  all  atoms  that  occurs  when  atom  a is  displaced,  and  the  term 
(overlap)  arises  primarily  from  the  change  in  the  "atomic  dipoles" 
from  the  electron  lone  pairs  as  the  a atom  moves,  as  well  as  from  the 
changes  in  atomic  dipole  contributions  from  the  bonding  electrons. 

The  "atomic  dipole"  was  defined  by  Coulson  (67)  as  the  contribution  to 
the  total  dipole  moment  due  to  the  displacement  of  the  centroid  of 
charge  due  to  the  electrons  from  the  centroid  of  positive  charge  due 
to  the  nucleus  of  atom  a.  By  examining  the  elements  of  all  these 
components  in  this  charge-charge  flux-overlap  (CCFO)  model  (66)  it 
is  possible  to  interpret  the  APTs  and  further  analyze  the  intensity 
parameters  of  a molecular  system.  Each  of  the  three  component  tensors 
also  has  invariant  properties.  The  mathematical  expressions  for  these 
quantities  were  developed  by  Chin  (28) . For  the  effective  charge 
of  the  charge  tensor  he  obtained 

X“(C)  = (1/3)^{p“(C)  V (3-15) 

where  Pa(C)i:L  is  a diagonal  term  of  the  charge  tensor  for  atom  a, 
given  by  (66) 

Pa(C) . . = Q = N - Z . (3-16) 

ii  a a a 

Here  Q is  the  Mulliken  net  atomic  charge,  N is  the  Mulliken 
a a 

population  (68),  and  Za  is  the  number  of  electrons  in  the  neutral  atom. 

The  effective  charge  x“(CF)  for  the  charge  flux  tensor  is  given  by 
Xa(CF)2  = Tr{pa(CF)[pa(CF^+}  (3-17) 

where  P(CF)  is  the  charge  flux  for  atom  a . 

We  are  interested  in  these  two  invariants  because  we  want  to 
obtain  the  charge  undeformability  values  for  the  APTs  of  uracil. 

This  parameter  was  defined  by  Chin  (28)  as 

charge  undeformability  = \a (C)  /xa (CF)  . (3-18) 
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It  was  suggested  by  Chin  (28)  that  this  ratio  is  a good  parameter  to 
measure  charge  flows. 

For  the  rest  of  the  invariant  properties  of  the  component  tensors 
the  reader  should  consult  Ref . 28 . 

Atomic  Polar  Tensors  for  Uracil 

The  equilibrium  geometry  of  uracil  used  throughout  all  this  study 
was  a slightly  modified  version  of  that  calculated  for  the  optimized 
geometry  by  Nishimura,  Tsuboi,  Kato,  and  Morokuma  (14)  using  the 
ST0-3G  basis  set.  The  C-H  and  N-H  bond  lengths  were  decreased  slightly 
from  their  predicted  values  to  bring  them  into  better  agreement  with 
experimental  values.  That  geometry  and  the  numbering  of  the  atoms 
used  in  these  calculations  are  given  in  Fig.  3-1.  In  addition,  Fig.  3-1 
indicates  the  orientation  of  the  principal  Cartesian  axes.  In  terms  of 
these  axes,  the  principal  Cartesian  coordinates  for  the  atomic  positions 
are  listed  in  Table  3-1.  The  calculated  values  for  the  principal 
moments  of  inertia  Ixx,  Iyy,  and  Izz  are  also  given  in  this  table.  The 
molecule  is  planar,  in  agreement  with  the  experimental  structure  found 
for  uracil  (69).  Therefore,  this  molecule  belongs  to  the  Cg  symmetry 
point  group. 

With  the  4-3 1G  basis  set  (70)  and  the  GAUSSIAN  76  system  of 
computer  programs  (71) , the  total  energy  calculated  for  uracil  with 
this  geometry  was  -411.85092268  Hartrees,  and  the  total  dipole  moment 
value  was  calculated  to  be  5.487  Debye,  with  components!  Px°  — 5.289  D, 
pyO  = _i.460  D,  and  pz°  = 0.000  D.  Experimentally,  Kulakowska,  Geller, 
Lesyng,  and  Wierzchowski  (72)  determined  that  the  total  dipole  moment 
for  uracil  in  dioxane  is  4.160  D,  within  experimental  uncertainty. 
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Figure  3-1.  Uracil  geometry  and  atom  numbering  as  used  in  the 
present  study.  The  molecule  is  oriented  in  its 
principal  Cartesian  coordinate  system. 
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TABLE  3-1 

PRINCIPAL  CARTESIAN  COORDINATES,  MOMENTS 
OF  INERTIA  AND  DIPOLE  MOMENT  COMPONENTS 
FOR  THE  ASSUMED  EQUILIBRIUM  GEOMETRY 
OF  URACIL 

Principal  Cartesian  Coordinates  (&) 


Atom 

X 

Y 

Z 

N1 

1.022080 

-1.201760 

0.000000 

c2 

-0.358920 

-1.203150 

0.000000 

n3 

-0.939870 

0.050810 

0.000000 

C4 

-0.292760 

1.270810 

0.000000 

c5 

1.144910 

1.186370 

0.000000 

c6 

1.723390 

-0.017890 

0.000000 

h7 

1.516120 

-2.087260 

0.000000 

°8 

-1.032580 

-2.215490 

0.000000 

«9 

-1.953500 

0.078090 

0.000000 

°10 

-0.959190 

2.298670 

0.000000 

H11 

1.722080 

2.115730 

0.000000 

h12 

2.790810 

-0.153530 

0.000000 

Principal  Moments  of  Inertia  (ui^) 
Ixx  = 245.860840 
Iyy  = 129.714401 
Izz  = 375.575195 

Dipole  Moment  Component 
px  = 5.289  D = 1.233  ei 
py  = -1.460  D = -0.304  el 
pz  = 0.000  D = 0.000  el 


aAtom  numbering  as  given  in  Fig.  3-1. 

k Experimental  dipole  moment  value  is  4.160  D;  Ref.  72. 
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The  APIs  for  uracil  were  calculated  by  Chin  according  to  procedures 
described  in  Chapter  2,  and  published  in  Ref.  17.  However,  a sign  mis- 
take was  made  while  writing-up  one  of  the  dipole  derivatives  in  the 
computer,  resulting  in  small  changes  in  some  of  the  APTs  reported  (17). 
The  implications  of  this  sign  change  will  be  discussed  in  Chapter  4, 
when  we  present  the  predicted  intensities  for  uracil.  The  calculated 
values  of  the  dipole  moment  derivatives  P^  are  displayed  in  Table  3-2. 
The  problem  just  mentioned  with  the  sign  was  in  the  x component  of 
the  i>2Q  coordinate.  Its  value  should  be  -0.0834  e.  Also,  the  values 
of  the  Pp  tensor  (Eq.  (2-37))  are  given  in  this  table. 

The  corrected  APTs  for  all  12  atoms  of  uracil  are  given  in 
Table  3-3  for  the  coordinate  system  shown  in  Fig.  3-1.  Note  that  the 
units  of  the  dipole  derivatives  in  the  APTs  are  in  multiples  of  the 
charge  of  the  electron  e,  on  atomic  units.  Conversions  from  other 
units  are  4.8  D A-1  = 1 e = 1.602  x 10-19  C.  The  APT  invariants  xa» 
p , and  were  calculated  using  Eqs . (3-6),  (3-7)  and  (3-8),  and  are 
also  included  in  Table  3-3.  The  xa  values  have  been  partitioned  into 
xa'  (for  in-plane)  and  x^"  (for  out-of-plane).  The  importance  of  these 
Xa  values  will  be  made  clear  when  we  use  them  in  the  next  section  with 
the  total  absolute  intensity  sum  for  uracil. 


TABLE  3- 
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CORRECTED  ATOMIC  POLAR  TENSORS  FOR  URACIL 
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Total  Absolute  Infrared  Intensity  Prediction  for  Uracil 

Since  now  the  APTs  and  dipole  moment  components  (Px°  * Py° » and 

p °)  are  given  in  Table  3-3,  and  the  principal  moments  of  inertia  in 
z 

Table  3-1,  the  intensity  sum  rule  (Eq.  (3-13))  can  be  applied  to  see 

what  is  predicted  for  the  ratio  of  the  total  intensity  sum  of  all 

the  in-plane  vibrations  (A')  to  the  total  sum  for  the  out-of-plane 

vibrations  (A")  of  uracil.  This  will  be  some  measure  of  the  accuracy 

of  the  prediction,  particularly  for  the  out-of-plane  modes.  For  A' 

vibrations  (active  in  the  xy  plane,  according  to  Fig.  3-1)  the  square 

of  the  effective  charge  (xa')2  for  atom  a will  be  given  by  (64) 

(Xaa'>2  = (l/3)[(3px/3Xa)2  + (3px/3Ya)2  + (3py/3Xa)2  + (3py/3Ya)2]  . 

(3-19) 

Also,  for  the  A"  vibrations  (active  in  the  xz  or  yz  planes  of 

Fig.  3-1)  the  square  of  the  effective  charge  (xa")  for  each  a atom  is 

(Xaa")2  = (l/3)(3pz/9Za)2  . (3-20) 

The  results  of  the  calculated  values  for  xa2>  (Xaa  )2>  and  (Xaa  )2  are 
displayed  in  Table  3-4. 

The  rotational  correction  expressions  ^a  and  0a  for  A'  and  A" 
vibrations,  respectively,  are  obtained  by  partitioning  the  result 
of  the  solution  of  Eq.  (3-13)  with  respect  to  the  xy  (for  A'),  and 
yz  or  xz  (for  A")  planes.  This  was  purposely  done  in  Eq.  (3-12) 
where  the  first  bracket  corresponds  to  0a 

aa'  - (p2°)2/Ixx  + <PZ°)2/Iyy  + t(Px°>2  + (Py°)2I/Izz  <3-21) 

ii 

and  the  second  bracket  gives  fta 

fia"  = (Py°)2/Ixx  + (P  X°)2/Iyy  • (3~22) 

These  rotational  correction  expressions  were  calculated  for  uracil, 
and  the  values  are  fla  =3.39  km  mol-^  and  fia  =9.49  km  mol 


TABLE  3-4 

INTENSITY  SUM  RULE  APPLIED  TO  URACIL 
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Predicted  intensity  : 2092.30  656.35  2748.65 

1 1 
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Accordingly,  the  sum  in  the  right-hand-side  of  Eq.  (3-13)  predicts 

for  A'  vibrational  modes  a total  intensity  of  2092.30  km  mol  * , while 

a total  sum  of  656.35  km  mol  ^ is  predicted  for  A"  vibrational  modes. 

Consequently,  the  ratio  between  the  predicted  total  intensity  sums  for 

in-plane  modes  to  out-of-planes  is:  2a  /Sa  = 3.19.  The  experi- 

a'  s a s 

mental  value  for  this  ratio  was  found  to  be  about  3.70,  from  measu- 
rement of  intensities  of  uracil  in  an  argon  matrix  (16,17). 


CHAPTER  4 

RESULTS  OF  EXPERIMENTAL  AND  CALCULATED 
SPECTRA  FOR  URACIL 

The  principal  considerations  in  this  chapter  are  the  normal  coordi- 
nate analysis  and  the  prediction  of  frequencies,  as  well  as  absolute 
infrared  intensities,  for  uracil.  We  shall  also  review  the  assignment 
of  the  infrared  spectrum  of  isolated  uracil  and  contrast  it  with  the 
spectrum  in  the  solid  phase,  in  solution,  etc.  In  order  to  carry  out 
these  objectives,  comparisons  are  made  between  our  force  constants  and 
assignments  for  uracil  with  those  reported  in  the  literature.  The 
importance  of  these  comparisons  is  to  illustrate  and  establish  the 
basis  for  future  studies  of  larger  molecules. 

Normal  Coordinate  Analysis  for  Uracil 

To  carry  out  the  normal  coordinate  analysis  (NCA)  of  uracil  we 
used  the  procedure  described  in  Chapter  2.  Since  uracil  has  12  atoms 
and  is  planar,  it  belongs  to  the  Cs  symmetry  group;  there  will  be  21 
in-plane  vibrations  of  A'  species,  and  9 out-of-plane  vibrations  of 
A"  species,  for  a total  of  30  vibrational  fundamentals. 

In  the  present  study  our  efforts  have  been  directed  to  the  21 
in-plane  vibrations  of  uracil,  for  the  following  reasons.  First, 
over  the  past  few  years  NCAs  and  assignments  of  the  infrared  spectrum 
of  uracil  have  been  carried  out  mainly  for  these  vibrations  (12—15), 
using  force  constants  obtained  employing  semiempirical  methods  (as  will 
be  discussed  in  more  detail  in  the  next  chapter) . Therefore,  we  can 
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make  comparisons  of  our  results  using  ab  initio  force  constants  with 
those  obtained  in  the  past,  and  evaluate  the  reliability  of  these  force 
constants.  Second,  understanding  uracil  in-plane  vibrations  should 
help  us  to  the  proper  interpretation  of  its  complex  spectrum,  as  well 
as  of  the  spectra  of  related  molecules,  such  as  the  methylated  uracils . 

Of  course,  the  out-of-plane  vibrations  are  of  interest  also,  but  more 
difficult  to  study  because  coordinates  for  the  out-of-plane  and  tor- 
sional modes  have  to  be  defined,  and  the  corresponding  force  constants 
are  less  well  known.  Besides,  the  quantum  mechanical  calculations 
have  not  yet  been  made  to  obtain  the  ab  initio  force  field  for  these 
vibrations . 

In  this  study  the  force  field  given  by  Nishimura,  Tsuboi,  Kato, 
and  Morokuma  (14)  was  used.  In  their  _ab  initio  calculation  they 
employed  the  ST0-3G  basis  set  to  obtain  energy  gradients  in  the  Cartesian 
coordinate  system  which  were  then  numerically  differentiated  in  order  to 
produce  the  force  constants  (73) . It  is  well  known  that  the  force  con- 
stants calculated  using  a 4-31G  basis  set  will  give  a better  reproduction 
of  the  observed  vibrational  frequencies  than  will  those  calculated  with 
the  ST0-3G  set.  Therefore,  Nishimura  et  al.  (14)  proposed  that  the 
"ST0-3G  force  constants"  be  multiplied  by  scaling  factors  chosen  on  the 
basis  of  their  experience  with  other  molecules  to  change  these  values  to 
values  for  the  force  constants  that  would  have  been  calculated  using  a 
4-31G  basis.  In  mathematical  terms  this  would  be  equivalent  to 

F (4-31G)  = c. .c. .F  (ST0-3G)  (4-1) 

ij  ii  3 J 13 

where  F (4-31G)  is  a 4-31G  force  constant,  F..(ST0-3G)  the  same  element 

ij 

obtained  using  the  ab  initio  ST0-3G  basis  set,  and  c^  and  c^ ^ are  the 
scaling  factors  corresponding  to  the  i and  j diagonal  force  constants. 
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If  the  element  is  diagonal  then 

F..(4-31G)  = CiiFi;L(ST0-3G)  . (4-2) 

For  example,  they  proposed  for  either  a C=0  or  C=C  stretch  that  the 
scaling  factor  is  0.82;  for  either  a C=0  or  C=C  bend  the  scaling  factor 
is  1.08,  etc.  (see  Table  I in  Ref.  14). 

The  internal  coordinates  for  uracil  are  defined  in  Fig.  4-1  and  the 
in-plane  symmetry  coordinates  are  given  in  Table  4-1.  Here,  the  redun- 
dancies of  the  angular  internal  coordinates  of  the  ring  system  have  been 
removed.  As  seen  in  Fig.  4-1,  there  are  30  internal  coordinates  defined 
to  describe  the  21  in-plane  degrees  of  freedom.  Hence,  the  number  of 
redundancies  is  9.  They  may  be  classified  as  follows:  six  associated 
with  the  angles  common  to  each  ring  atom,  for  example  A S 3 + Ag^g  + 

^®25  = 0 (see  Fig.  4-1);  one  cyclic  redundancy  corresponding  to  the 
constraint  that  all  internal  angles  of  the  ring  cannot  be  increased 
simultaneously,  Agj^  + Ag^  + Ag-^  + Agjg  + Ag^y+  Ag^g  = 0;  and  two 
other  cyclic  constraints  which  express  the  condition  that  the  ring  is 
closed  so  that  linear  combinations  of  the  stretching  and  bending  internal 
coordinates  are  constrained  to  zero  (6,21).  The  removal  of  these  nine 
redundancies  makes  the  calculation  easier,  since  they  can  just  be 
eliminated  after  constructing  the  remaining  coordinates  to  be  orthogonal 
to  the  redundancy  condition.  For  this  reason  the  bending  coordinates  in 
Table  4-1  appear  as  linear  combinations  (73)  . Based  on  the  symmetry 
coordinates  of  Table  4-1,  the  in-plane  ab  initio  force  constants  for 
uracil,  scaled  to  give  "4-31G"  values,  are  given  in  Table  4-2.  Using 
these  force  constant  values  we  carried  out  the  NCA,  setting  up  the  B and 
G matrices  (see  Chapter  2) . The  eigenvalues  were  obtained  and  the  fre- 
quencies calculated  by  solving  Eq.  (2-21).  We  found  that  the  calculated 
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Hl2 


Figure  4-1.  In-plane  internal  coordinates  of  uracil. 
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TABLE  4-1 

IN-PLANE  SYMMETRY  COORDINATES  FOR  URACIL 


S1  -4rl 

vN^C2  stretching 

S2  ' 4r2 

vC2N3  stretching 

S3  -Ar3 

vN^C^  stretching 

S4  4r4 

vC^C,.  stretching 

S5  -4r5 

vC^Cg  stretching 

S6  '4r6 

vC,N.  stretching 
0 I 

S7  = 4r7 

vN^H^  stretching 

S8  ' 4r8 

vC20g  stretching 

S9  Ar9 

vNgHg  stretching 

S10  = 4r10 

vC4°10  stretcllin§ 

Sn  = 4rn 

vC^H^  stretching 

S12  " 4r12 

vC6H12  stretchin8 

S13  ‘ (4613-4614+4815-4ei6+4S17-4S18)/^ 

6Ring  deformation 

I 

Su  = ( AS13-Ae14-A815+  A816-A617-A813)/^ 

6Ring  deformation 

II 

S15  ' <4814-4S15+4S17-4S18)/2 

<$Ring  deformation 

III 

S16  * WB19-AS25)/^ 

6N^H  bending 

S17  ’ <4820-4e26,/^ 

5C20  bending 

S18  = (A621-AB27)/^2 

6NgH  bending 

S19  = (A322-AB28)/^2 

6C^0  bending 

S20  ‘ (4823-4629)/’^ 

6C^H  bending 

S21  ■ (4824-4830>/^ 

6C,H  bending 
0 

Note:  The  atomic  numbering  is  shown  in  Fig.  4-1,  which  also  gives 
the  internal  coordinate  definitions.  Since  the  only  element 
of  symmetry  (besides  E)  is  each  internal  stretching 

coordinate  is  a symmetry  coordinate.  But  to  represent 
bending  motions,  with  internal  coordinates,  appropriate 
linear  combinations  must  be  made. 


TABLE  4- 
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'The  S symmetry  coordinates  correspond  to  those  given  in  Table  4-1. 

'Units:  md  A”*  for  stretch-stretch  elements,  md  A for  bend-bend  elements,  and  md  for  bend-stretch  elements. 

These  force  constant  values  are  obtained  using  the  scaling  factors  given  by  Nishimura,  Tsuboi,  Kato,  and  Morokuma  (14)  to  the  ST0-3G 
force  field  that  they  obtained. 
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frequencies  should  be  further  scaled  by  0.90  to  bring  them  into  better 
agreement  with  experimental  data.  This  scaling  is  necessary  for 
frequencies  calculated  using  the  4-31G  force  field.  For  example,  a 
similar  conclusion  was  reached  by  Wiberg,  Walters,  Wong,  and  Colson  (74), 
in  their  study  of  frequencies  and  intensities  of  pyridine.  The 
resulting  adjusted  frequencies  for  the  in-plane  normal  modes  of  uracil 
are  given  and  compared  with  the  observed  values  in  Table  4-3,  which  also 
includes  contributions  greater  than  5 % to  the  potential  energy  distri- 
butions (PEDs,  see  Eq.  (2-29))  for  each  normal  mode.  The  frequencies 
that  we  calculated  agreed  with  predictions  by  Nishimura,  Tsuboi,  Kato, 
and  Morokuma  (14)  for  uracil  within  3 %.  This  slight  discrepancy  is 
due  to  the  small  differences  in  the  G matrix,  arising  from  the  slight 
differences  in  our  assumed  geometry  as  compared  to  theirs.  The  results 
of  our  study  have  been  reported  and  discussed  by  Chin,  Scott,  Szczepaniak, 
and  Person  (17).  Our  assignments  of  the  uracil  spectrum  will  be 
discussed  in  the  next  sections. 

Comparison  of  Results 

In  this  section  we  want  to  establish  a comparison  between  the 
ab  initio  force  constants  for  uracil  and  the  semiempirical  force  con- 
stants reported  in  the  literature.  We  are  also  interested  in  com- 
paring the  calculated  frequencies  using  both  ab  initio  and  semiempirical 
constants,  as  well  as  the  assignment  descriptions  or  PEDs.  Our  motive 
for  these  comparisons  is  twofold:  first,  this  will  give  us  some  insight 
into  the  reliability  of  ab  initio  force  constants,  and  second  we  will 
verify  that  our  calculation  and  assignment  of  spectra  gives  expected 
results.  This  understanding  should  prove  of  some  value  when  we  study 


larger  molecules. 


TABLE  4-3 

RESULTS  FOR  URACIL  IN-PLANE  VIBRATIONS 
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1389  21  1382  4 6CgH(29),  6C5H(15),  6NjH(14) , 6N3H(14),  vC40(8),  vC20(7) 

1219  4 1254  18  vN3C4(27) ,6N3H(22) , 6C5H(14),  vC2N3(13),  vC4C5(ll),  vNjC2(10) , 6C40(10) 

1186  109  1175  44  6CgH(40),  6C5H(22),  6KjH(15) , vCgNj(9) 
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Force  constants  in  uracil  studies  compared 

In  general,  force  constants  obtained  by  ab  initio  quantum  mechanical 
methods  show  systematic  deviations  from  those  considered  the  most  relia- 
ble ones  determined  semiempirically  from  a large  number  of  independent 
experimental  data  by  spectral  fit  refinements,  as  will  be  discussed  in 
the  next  chapter.  We  are  mainly  concerned  with  the  patterns  followed  by 
the  force  constants  with  the  use  of  each  approach  (ab  initio  and  semi- 
empirical)  . Some  illustrative  comparisons  of  force  constants  for  uracil 
are  given  in  Table  4-4.  The  first  column  shows  values  for  some  of  the 
diagonal  force  constants  obtained  by  Nishimura,  Hirakawa  and  Tsuboi  (6, 
to  which  we  will  refer  as  NHT) . The  treatment  they  followed  involved  a 
spectral  fit  refinement  using  the  experimental  frequencies  obtained  for 
uracil  and  three  deutero-derivatives.  They  attempted  to  reproduce  the 
observed  frequencies  for  these  molecules  carrying  out  a NCA  with  force 
constants  obtained  by  means  of  an  iterative  procedure.  Column  2 refers 
to  force  constants  from  a simple  valence  force  field  calculated  by  Susi 
and  Ard  (12,  SA) . They  used  frequencies  from  four  isotopic  uracil 
species,  for  a total  of  83  experimental  frequencies,  in  their  refinement 
to  obtain  23  force  constants.  Column  3 corresponds  to  the  study  by 
Bandekar  and  Zundel  (3,  BZ)  on  crystalline  uracil.  Initially,  they  con- 
structed a force  field  for  uracil  transferring  force  constants  from 
small  model  compounds,  such  as  N-methylacetamide  and  acetone.  The  final 
set  of  force  constants  involved  the  frequencies  of  all  four  uracil 
molecules  for  the  complete  unit  cell  (48  atoms) . Force  constants  of  the 
study  by  Bowman  and  Spiro  (15,  BS)  are  given  in  column  4.  They  calcu- 
lated a force  field  for  uracil  using  the  semiempirical  molecular 
orbital  calculation,  MNDO,  and  a method  of  molecular  orbital 


TABLE  4-4 

COMPARISON  OF  DIAGONAL  FORCE  CONSTANTS  FOR  URACIL  REPORTED  IN  THE  LITERATURE 
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constraint  on  the  interaction  coordinates.  Then,  they  evaluated  a final 
set  of  force  constants,  starting  from  their  initial  MNDO  values,  via  a 
nonlinear  least-squares  fit  of  the  diagonal  elements  to  observed  frequen- 
cies. Also  in  Table  4-4,  column  5 presents  some  force  constants  obtained 
using  the  semi- empirical  CNDO/2  method  by  Harsanyi  and  Czaszar  (75,  HC) . 
They  made  use  of  correcting  factors  (obtained  from  their  studies  of 
maleimide)  for  these  CNDO/2  values  to  reproduce  the  gas  phase  infrared 
spectrum.  In  column  6 we  show  some  values  obtained  at  the  ST0-3G  level 
and  corrected  to  "4-3 1G  level",  in  the  calculation  by  Nishimura,  Tsuboi, 
Kato,  and  Morokuma  (14,  NTKM) . These  are  the  same  force  constants  that 
we  are  using,  and  described  previously,  in  the  present  study. 

An  inspection  of  Table  4-4  shows  that  the  empirically  obtained 
force  constant  values  (NHT,  SA,  and  BZ)  are  very  similar.  The  same 
does  not  hold  true  for  those  found  using  different  theoretical  methods 
(BS,  CH,  and  NTKM).  Nevertheless,  the  pattern  followed  by  ab  initio 
force  constants  compares  very  well  with  that  followed  by  the  empirical 
ones . 

Let  us  examine  some  particular  cases  in  the  comparison  of  force 
constants  given  in  Table  4-4.  The  C2=0  stretching  force  constant  values 
are:  11.0  (NHT),  11.0  (SA) , and  9.62  (BZ)  mdyn  for  the  empirical 
cases,  and  9.601  (BS) , 10.131  (HC) , and  13.298  (NTKM)  mdyn  ft-1  for  the 
calculated  values.  The  NTKM  value  of  13.298  mdyn  X-1  is  much  higher 
than  for  all  other  instances.  Since  we  must  scale  the  calculated 
frequencies  by  a factor  of  0.90  to  get  agreement  with  observed  data, 
all  4-31G  force  constants  (NTKM)  must  be  scaled  by  0.81;  the  results  of 
this  scaling  are  presented  in  column  7 of  Table  4-4.  After  scaling. 
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the  value  of  13.298  mdyn  & ^ for  the  stretching  force  constant  by 

NTKM  is  reduced  to  10.771  mdyn  X which  agrees  very  well  with  the 
other  ones  listed  above. 

For  the  N^H  bend,  the  force  constant  values  in  Table  4-4  are: 
empirically,  0.510  (NHT) , 0.417  (SA) , and  0.510  (BZ) ; calculated, 

0.641  (HC)  and  0.619  (NTKM).  Again,  those  values  calculated  quantum 
mechanically  are  higher  than  the  empirical  ones.  The  value  from  NTKM 
after  scaling  by  0.81  of  0.501  mdyn  2.  is  similar  to  those  obtained 
empirically. 

It  is  not  possible  to  make  similar  comparisons  for  many  interaction 
constants  because  what  was  determined  in  the  past  empirically  for  a 
molecule  were  simple  valence  force  fields  containing  only  diagonal  terms. 
Nearly  all  other  interaction  constants  had  to  be  assumed  to  be  zero. 

The  interaction  force  constants  for  the  C^=0  stretch/C^=C^  stretch  were 
however  obtained  in  previous  studies  and  are  given  in  the  last  line  of 
Table  4-4.  All  we  can  say  is  that  its  value  was  expected  to  be  small, 
and  probably  negative,  as  given  by  SA  (-0.225  mdyn  & ^),  HC  (-0.090 
mdyn  &-1),  and  NTKM  (-0.128  or  -0.104  mdyn  ft-1  when  scaled  by  0.81). 

The  evidence  presented  in  Table  4-4  supports  the  idea  that  ab  initio 
force  constants  are  reliable  enough  to  be  used  in  NCA  calculations.  We 
strongly  recommend  their  use  when  proper  scaling  factors  are  applied  to 
them,  since  our  results  (17-19)  clearly  show  that  these  force  constants 
reproduce  the  experimental  spectrum,  within  + 20  cm  in  the  frequencies. 
This  value  was  obtained  through  a least  squares  procedure  of  the 
experimental  and  calculated  frequency  data  for  uracil  (17)  and 
methylated  uracils  (18,19). 
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Comparison  of  Infrared  Spectral  Assignments  for  Uracil 

Now  we  may  compare  the  sets  of  frequencies  calculated  by  NCA  calcu- 
lations by  the  different  workers  mentioned  above.  These  sets  of  values 
are  shown  in  Table  4-5.  First,  it  is  worth  noting  just  what  experimental 
data  were  used  in  the  fitting  of  frequencies  in  each  study.  It  is 
briefly  summarized  in  the  following:  NHT  (6),  resonance  Raman  spectra  of 
neutral  aqueous  solutions  of  uracil  plus  results  from  SA;  SA  (12),  ordi- 
nary laser  Raman  and  infrared  spectra  for  polycrystalline  uracil;  BZ  (3), 
also  infrared  and  Raman  spectra  for  solid  uracil;  BS  (15)  used  the 
experimental  frequencies  of  SA;  and  HC  (75)  tried  to  fit  data  from  the 
gas  phase  infrared  spectrum  of  uracil.  Actually,  NTKM  (14)  did  not  try 
to  fit  the  calculated  frequencies  to  the  observed  values,  but  instead  to 
see  how  close  the  predicted  frequencies  using  ab  initio  force  constants 
are  to  the  observed  ones.  In  the  present  study  we  used  predicted  fre- 
quencies and  infrared  intensities,  as  well  as  PEDs,  to  decide  which 
experimental  band  in  the  matrix-isolated  spectrum  of  uracil  was  to  be 
assigned  to  which  calculated  mode. 

A close  examination  to  Tables  4-3  and  4-5  shows  a surprising  simi- 
larity in  the  frequency  assignment  descriptions  in  terms  of  the  PEDs. 

For  example,  except  for  the  BZ  assignment  of  the  two  N-H  stretches,  in 
which  they  calculated  about  half-and-half  mixing  of  these  modes,  the 
rest  of  the  workers  agree  that  the  two  N-H  stretches  are  local  modes. 

The  spectral  region  in  which  the  C=0  stretches  and  C =C,  stretch  are 

j o 

assigned  is  rather  interesting.  Here,  the  general  consensus  is  that 
the  C2=0  stretching  frequency  is  higher  than  the  corresponding  value 
for  the  C^=0  stretch.  The  calculated  listed  values  in  Table  4-5 
going  from  left  to  right  are  for  C2O  : 1697,  1684,  1721,  1725,  1765 


TABLE  4-5 

COMPARISON  OF  FREQUENCIES  FOR  URACIL  FROM  DIFFERENT  STUDIES 
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and  1775  cm  and  for  C,0  are  : 1648,  1653,  1673,  1669,  1739,  and 

4 

1748  cm  . It  should  be  remembered  that  the  differences  in  these  calcu- 
lated values  are  mainly  due  to  attempts  to  fit  spectral  data  obtained  in 
different  phases  (aqueous,  solid,  etc.),  where  hydrogen  bonding  lowers 
this  mode.  Anyway,  the  PEDs  reported  by  NHT,  SA,  and  HC  compare  very 
well  for  these  modes  with  those  of  our  study  (Table  4-3).  The  value 
predicted  for  the  C^=Cg  frequency  and  PEDs  are  almost  the  same  in  all 
the  studies,  except  BS,  who  do  not  predict  it  to  be  important  in  any 
mode . 

The  remaining  assignments  for  the  uracil  spectrum  are  similar  in 
all  the  seven  studies,  as  shown  in  Tables  4-3  and  4-5.  However,  many 
of  the  modes  show  extensive  mixing  which  is  sensitive  to  the  values  of 
the  force  constants.  Everything  considered,  the  agreement  between 
these  different  treatments  is  remarkable,  but  the  differences  are 
important . 

Predicted  Absolute  Infrared  Intensities  for  Uracils 

The  predicted  infrared  intensities  for  uracil  and  seven  of  its 
deuterated  derivatives  were  reported  by  Chin,  Scott,  Szczepaniak,  and 
Person  (17).  Since  we  had  to  correct  the  APTs  for  uracil  (as  explain- 
ed in  Chapter  3)  it  is  now  appropriate  to  describe  the  effect  of  such 
correction  on  the  predicted  intensities.  First,  the  only  change  occurs 
in  the  intensities  of  the  following  modes  (see  Table  I of  Ref.  17): 
v^(NjH  stretch),  V2(N^H  stretch),  v^(mainly  C2=0  stretch),  and 
(mainly  C^=0  stretch) . For  the  PED  description  of  these  modes  see 
Table  4-3. 

In  Table  4-6  are  given  the  wavenumbers  and  calculated  intensities 
for  the  four  vibrational  modes  just  mentioned.  Interestingly,  the 


TABLE  4-6 

CORRECTED  INTENSITIES  FOR  MODES  v, , v0,  v,,  AND  v,  OF  URACILS 
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mode  for  uracil  now  has  an  intensity  of  108  km  mol-* , and  V2  a value 
of  97  km  mol-* . These  results  now  agree  better  with  the  experimentally 
observed  intensities.  It  should  be  mentioned  that  previously  (17)  one 
of  the  dramatic  differences  between  our  calculated  spectrum  and  the 
experimental  matrix  spectrum  was  the  predicted  "reversal"  in  this 
pattern  of  intensities.  Another  important  result  is  evident  from 
Table  4-6  concerning  the  predicted  intensities  of  and  . For  uracil 
these  values  are  respectively  993  and  458  km  mol-*,  so  that  the  calcu- 
lated ratio  of  the  intensities  A^/Ag  is  2.17,  while  previously  it  was 
1.64  (17).  The  experimental  intensity  ratio  for  the  same  modes  is  about 
2.34,  taking  the  ratio  of  the  values  of  the  relative  integrated  absorb- 
ances for  these  modes,  680/291  respectively,  reported  by  Chin  et  al. 

(17)  for  uracil. 

The  right  pattern  is  also  predicted  for  the  intensities  of  the  N-H 
and  C=0  stretching  modes  of  N-^  jN^-dideuterouracil,  as  shown  in  Table 
4-6  (compare  with  Table  II  of  Ref.  17).  The  N^D  stretching  mode  has  an 
experimental  relative  intensity  of  112,  and  the  corresponding  value  for 
the  N^D  stretch  is  73.  The  predicted  respective  values  given  in  Table 
4-6  are  81  and  73  km  mol-* . The  ratio  of  the  experimental  relative 
integrated  absorbance  of  the  C2=0  stretching  mode  to  that  of  the  C2=0 
stretching  mode  is  391/612,  or  0.64,  while  the  calculated  corresponding 
value  from  Table  4-6  is  479/930,  or  0.52. 

Our  assignment  (17)  of  the  vibrational  spectrum  of  uracil  does  not 
need  to  be  changed  on  the  basis  of  the  new  values  of  predicted  intensi- 
ties. The  new  corrected  results  for  the  predicted  intensities  give  us 
more  confidence  in  the  validity  of  the  calculated  APTs  for  uracil. 
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However,  it  seems  that  one  must  be  extremely  careful  in  carrying  out 
calculations  of  intensities  because  the  simplest  mistake  can  have 
important  implications. 

Simulation  of  Infrared  Spectra 

In  order  to  have  a visual  comparison  between  the  observed  spectrum 
and  the  calculated  frequencies  and  absolute  intensities  for  uracil,  the 
calculated  spectrum  was  simulated  and  compared  with  a corresponding 
simulation  of  the  matrix  spectrum.  For  the  simulation  of  the  calcu- 
lated spectrum,  all  bands  were  assumed  to  have  the  same  constant  half- 
width (5  cm  ^ , estimated  from  the  experimental  matrix  spectrum)  and  a 
peak  height  chosen  to  give  an  area  proportional  to  the  calculated 
intensity  for  each  band  from  each  normal  mode  at  the  calculated  fre- 
quency. The  intensity  proportionality  factor  was  chosen  empirically 
so  that  the  sum  of  all  the  intensities  for  the  in-plane  modes  in  the 
experimental  spectrum  is  equal  to  the  sum  of  intensities  from  the 
ab  initio  calculation  of  APTs,  using  the  intensity  sum  rule  to  obtain 
the  total  predicted  absolute  intensity,  as  described  in  the  last 
section  of  Chapter  2.  There  we  found  that  the  total  intensity  pre- 
dicted for  the  entire  spectrum  of  uracil  is  2749  km  mol  ^ (Table  3-4) . 

The  simulation  for  the  ab  initio  predicted  spectrum  for  the  21 
in-plane  normal  modes  of  uracil  is  shown  at  the  top  part  of  Fig.  4-2, 
along  with  the  simulated  experimental  spectrum  of  matrix-isolated 
uracil  (in  argon)  shown  at  the  bottom  part  of  the  figure,  which  includes 
out-of-plane  vibrations.  These  simulations,  and  a comparison  with  the 
spectrum  predicted  by  Harsanyi  and  Czaszar  (75)  for  uracil  using  the 
CNDO/2  method  have  been  discussed  by  Szczepaniak,  Szczesniak,  Nowak, 
Scott,  Chin,  and  Person  (18). 
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Figure  4-2.  Simulated  infrared  spectra  for  uracil,  a)  ab  initio 
predicted  spectrum;  b)  Matrix  isolation  infrared 
spectrum  in  argon  matrix.  (Redrawn  and  corrected  from 
Refs.  17  and  18). 
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Discussion  of  Results 

Examining  Table  4-3  and  Fig.  4-2  we  can  see  that  the  predicted  fre- 
quency and  absolute  intensity  pattern  agrees  closely  with  the  experi- 
mental matrix  spectrum  for  uracil,  within  the  typical  "factor-of-two" . 
This  excellent  agreement  leads  us  to  believe  that  the  ab  initio  force 
field  at  the  "4-3 1G  level"  that  was  used  in  the  NCA  is  a good  approxi- 
mation for  the  in-plane  modes  of  uracil. 

As  it  was  mentioned  in  the  previous  section,  Harsanyi  and  Czaszar 
(75)  studied  both  the  in-plane  and  out-of-plane  vibrations  of  uracil 
by  carrying  out  a NCA  with  a force  field  obtained  at  the  CNDO/2  level. 
This  force  field  was  then  scaled  by  correcting  factors  in  order  to 
yield  calculated  frequencies  in  good  agreement  with  the  experimental 
gas  phase  infrared  spectrum.  These  scaling  factors  were  taken  from 
studies  of  maleimide  (a  five  atom  molecule  with  some  similarity  to 
uracil),  and  then  used  with  the  rough  CNDO/2  force  field  calculated 
for  uracil.  Some  of  the  force  constant  values  so  obtained  were  given 
in  column  5 of  Table  4-4.  The  frequencies  calculated  with  this  force 
field  are  shown  in  column  5 of  Table  4-5. 

4 

We  have  compared  Harsanyi  and  Czaszar  (75)  results  with  our  cal- 
culation and  with  the  experimental  spectrum  (18).  We  found  that  the 
most  drastic  discrepancy  between  their  prediction  and  the  experimental 
spectrum  exists  in  the  NH  and  CH  stretching  regions  (see  Fig.  1A  of 
Ref.  18).  The  CNDO/2  calculation  predicts  much  stronger  intensities 
for  the  CH  stretches  than  for  the  NH  stretches.  This  result  contrasts 
dramatically  with  the  experimental  spectrum,  where  the  NH  stretches 
are  quite  intense  but  the  CH  stretches  are  hardly  visible.  Other  than 
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that,  the  scaled  CNDO/2  prediction  of  the  spectrum  agrees  surprisingly 
well  with  the  experimental  spectrum  and  with  the  more  complicated 
ST0-3G  (and  4-3 1G)  ab  initio  calculation. 

For  our  assignment  of  the  uracil  infrared  spectrum  we  have  relied 
heavily  on  the  experimental  spectrum  of  the  matrix  isolated  sample  ob- 
tained by  Szczesniak,  Nowak,  Rostkowska,  Szczepaniak,  Person,  and 
Shugar  (16)  (see  Chapter  2).  Very  recently,  Maltese,  Passerini, 
Nunziante-Cesaro,  Dobos,  and  Harsanyi  (2)  and  Barnes,  Stuckey,  and 
Le  Gall  (5)  have  also  reported  spectral  data  for  matrix  isolated 
uracils  from  studies  in  argon  matrices.  All  three  sets  of  data  for  ura- 
cil from  studies  in  Ar  matrices  have  been  collected  in  Table  4-7,  together 
with  the  proposed  assignment  of  the  infrared  spectrum  from  each  study. 

In  order  to  discuss  and  compare  the  assignment  for  the  uracil 
infrared  spectrum  we  have  partitioned  the  spectrum  in  three  regions: 
the  4000  - 2000  cm  ^ region,  the  2000  - 1600  cm  ^ region,  and  the 
1600  - 400  cm  ^ region. 

Spectral  Region  Between  4000  - 2000  cm  ^ 

The  highest  absorptions  in  this  region  are  assigned  to  the  two  NH 
stretches.  All  three  studies  presented  in  Table  4-7  agree  in  the 
assignment  of  the  3485  cm  ^ to  N^H  and  3435  cm  * to  N^H.  Slight 
discrepancies  in  frequency  suggests  some  calibration  errors.  There- 
fore, there  is  consensus  for  the  NH  stretching  region  assignment  by 
the  three  works. 

The  predicted  C^H  and  C^H  stretching  frequencies  of  3070  cm  * 

(11  km  mol  and  3031  cm  ^ (9  km  mol  for  uracil  (see  Tables  4-3 
and  4-7)  were  assigned  without  much  basis  (16)  to  very  weak  absorption 
bands  observed  at  3130  cm  ^ (Rel.  Int.  of  4)  and  2970  cm  ^ (Rel.  Int. 
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^The  PEDs  are  calculated  using  Eq.  (2-29).  See  Table  4-3  for  the  results  of  uracil. 
^Frequency  observed  in  the  Raman  spectrum  only. 
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of  8).  Barnes,  Stuckey,  and  Le  Gall  (5)  observed  a medium  Raman  band 
at  3084  cm  ^ , which  they  assign  to  the  C..H  stretching.  We  believe  this 
assignment  is  probably  correct.  The  C^H  stretch  was  not  observed  by 
Barnes  et  al.  (5)  for  uracil  in  the  argon  matrix,  but  for  deuterated 
uracils  (C^,C^-D  and  N^,N^-D)  in  the  solid  phase  at  20  K they  observed 
respectively  a "shoulder"  at  2295  cm  * and  a weak  band  at  3000  cm  ^ , 
which  they  assign  to  this  mode.  Neither  one  of  the  CH  stretching 
frequencies  were  observed  in  the  infrared  studies. 

Spectral  Region  Between  2000  - 1600  cm  ^ 

According  to  the  prediction  (see  Tables  4-3  and  4-7)  we  expect  to 
find  three  more  fundamental  modes  of  vibration  for  uracil  in  this  re- 
gion, namely  the  C^O  and  C^O  carbonyl  stretches  and  the  ring 

stretch.  We  predicted  the  C^O  stretching  mode  at  1775  cm  * (993  km 

mol  1),  the  C.O  stretch  at  1760  cm  ^ (458  km  mol  and  the  C =C, 

4 5 o 

stretch  at  1671  cm  * (2  km  mol  ^) . The  experimental  spectrum  shows 
that  this  is  a very  complex  region.  It  has  been  suggested  by 
Szczesniak,  Nowak,  Rostkowska,  Szczepaniak,  Person,  and  Shugar  (16) 
that  the  most  likely  explanation  for  such  complexity  is  that  Fermi 
resonance  occurs  between  the  strong  carbonyl  fundamentals  and  overtone 
or  combination  bands  from  the  lower  frequency  modes. 

The  assignments  by  the  three  studies  in  Table  4-7  assign  the  02=0 
stretching  mode  at  a higher  frequency  than  the  C^=0  one.  The  following 
observed  frequencies  in  cm  * are  assigned  to  the  02=0  and  C^=0  stretches 
respectively:  Barnes  et  al.  (5),  1764  vs/1761  vs  and  1706  vs;  Maltese 
et  al.  (2),  1763  vs/1761  vs/1756  ms/1740  ms/1733  ms/1727  s and 
1706  s/1703  s;  and  present  study,  1774  m/1762  vs/1758  s/1733  s and 
1720  w/1707  s/1699  mw.  Our  prediction  of  intensities  indicates 
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that  A,.  is  about  twice  A^;  in  fact  this  is  the  pattern  observed  in  the 
matrix  infrared  spectrum.  There  is  also  quite  general  agreement  in  the 
studies  concerning  the  assignment  of  the  stretching  mode,  except 

that  our  calculation  of  PEDs  shows  strong  mixing  of  this  mode  with 
other  modes.  However,  the  assignments  by  Barnes  et  al.  (5)  and 
Maltese  et  al.  (2)  consider  this  mode  to  have  a more  localized  character 
The  observed  frequencies  in  cm  ^ assigned  to  this  mode  are:  Barnes  et  al 
(5),  1643  m;  Maltese  et  al.  (2),  1643  m;  and  the  present  study,  1644  mw. 
Our  intensity  calculation  for  this  mode  is  2 km  mol  while  its 
relative  experimental  value  is  33. 

Spectral  Region  Between  1600  - 350  cm  * 

The  remaining  14  in-plane  fundamentals  of  uracil  are  located  in 
this  region  (see  Tables  4-3  and  4-7) . An  inspection  of  Table  4-7 
shows  that  there  is  a close  agreement  between  the  observed  frequencies 
of  Barnes  et  al.  (5)  and  those  observed  by  Szczesniak  et  al.  (16)  for 
uracil  in  an  Ar  matrix.  However,  there  are  some  discrepancies  with 
the  frequencies  observed  by  Maltese  et  al.  (2),  especially  the  lower 
values.  In  particular,  it  looks  as  if  the  587  cm  ^ (vw)  value  of 
Maltese  et  al.  (2)  was  caused  by  some  impurity,  because  if  this  value 
were  eliminated  and  the  last  three  frequencies  placed  up  by  one  step, 
then  there  would  be  a good  agreement  between  all  three  studies. 


CHAPTER  5 

TRANSFER  OF  VIBRATIONAL  PARAMETERS 


In  this  chapter  we  want  to  discuss  how  feasible  or  realistic  it  is 
to  predict  frequencies  and  infrared  intensities  of  different  molecules 
making  use  of  some  at  priori  knowledge  of  the  force  field  and  APTs  of 
structurally  and  chemically  related  molecules.  First  we  shall  consider 
the  problem  in  general.  Some  examples  will  be  given  concerning  the 
transfer  of  the  vibrational  parameters — force  constants  and  APTs.  For 
the  particular  case  of  uracil  a detailed  comparison  will  be  made  of 
these  two  sets  of  parameters  with  the  values  obtained  for  other  mole- 
cules. Evidence  supporting  the  hypothesis  that  the  transfer  of  force 
constants  and  APTs  is  a valid  procedure  will  lead  us  to  apply  it  to  the 
transfer  from  uracil  to  1-methyluracil,  3-methyluracil,  and 
1 , 3-dime thy luracil . 

Transfer  of  Force  Constants 

The  term  "force  constant"  is  used  specifically  for  the  harmonic 
(quadratic)  constants.  Except  for  the  simplest  molecules,  the  task  of 
determining  force  constants  from  experimental  data  is  a difficult 
problem,  due  to  the  large  number  of  unknown  constants  compared  to  the 
amount  of  available  data.  For  a secular  equation  of  order  n,  there 
are  (l/2)(n)(n  + 1)  independent  force  constants,  but  there  are  only  n 
fundamental  vibration  frequencies  (21,24-26,76).  Therefore,  it  is  in 
general  impossible  to  calculate  a unique  set  of  values  for  those  force 
constants  from  a set  of  frequencies  alone.  In  the  case  of  large 
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polyatomic  molecules,  Zerbi  (77)  has  pointed  out  that  since  no  unique 
choice  of  the  vibrational  potential  can  be  indicated,  only  chemical 
and  physical  common  sense  as  well  as  a critical  examination  of  the 
results  for  derived  physical  quantities  associated  with  vibrational 
motions  may  be  suggested  as  general  criteria  for  wise  handling  of  the 
force  constants. 

In  order  to  obtain  a physically  meaningful  set  of  force  constants, 

therefore,  one  should  employ  as  many  additional  data  for  these  physical 

quantities  associated  with  the  vibrational  motions  (such  as  isotope 

shifts  of  frequencies,  values  of  Coriolis  constants,  mean  amplitudes 

of  vibration,  etc.)  as  is  possible.  The  small  shift  in  vibrational 

frequencies  caused  by  isotopic  substitution  are  sensitive  functions  of 

the  force  field  (78)  since  the  potential  remains  the  same  from  one 

isotopic  species  to  another.  For  example,  the  study  of  Chalmers  and 

13 

McKean  (79)  on  CF  and  CF  demonstrates  that  quite  accurate  force 
4 4 

and  frequencies  are  obtained  by  the  use  of  isotopic  data. 

Another  valuable  source  of  information  in  the  force  constants 
determination  are  the  values  of  the  Coriolis  constants  (80)  which 
relate  the  mutual  interactions  of  the  normal  vibrations  through 
molecular  rotation.  Physically,  these  constants  measure  the  vibra- 
tional angular  momentum  generated  by  the  movement  of  the  atoms  in  the 
vibrations.  They  are  of  particular  value  in  determining  the  values  of 
the  interaction  force  constants  connecting  the  symmetry  coordinates 
used  (42,43).  However,  precise  evaluation  of  the  Coriolis  constants 
for  a molecule  is  very  difficult,  particularly  for  a large  molecule 


like  uracil. 
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Additional  data  that  may  be  used  to  evaluate  force  constants 
include  mean-square  amplitudes  of  vibration,  which  are  determined 
experimentally  as  part  of  the  structure  determination  by  the  gas-phase 
electron  diffraction  technique  (81).  Their  relationship  to  the 
molecular  force  field  arises  because  their  magnitudes  depend  on  the 
squares  of  the  relative  Cartesian  displacement  in  the  normal  vibrations 
of  different  pairs  of  atoms  (82) . 

Further  experimental  data  such  as  the  inertia  defect  parameters 
for  planar  molecules  (83),  the  £-type  doubling  constants  for  symmetric 
and  spherical  tops  (84),  etc.,  may  be  also  employed  when  attempting  to 
determine  experimental  molecular  force  fields. 

The  most  used  approach  toward  the  evaluation  of  the  force  con- 
stants for  large  molecules  involves  some  method  of  adjusting  the  force 
constants  so  that  the  calculation  reproduces  the  values  of  the  larger 
number  of  independent  experimental  parameters  as  closely  as  possible 
(26,85-95).  Generally,  the  method  of  least-squares  refinement  is 
used  for  this  purpose.  In  such  an  iterative  refinement  procedure 
one  modifies  a "guessed"  force  field  (obtained  by  transferring  force 
constants  from  similar  fragments  in  smaller  molecules)  to  fit  a gen- 
erally larger  number  of  experimental  frequencies  and  other  parameters. 
Then  the  calculated  eigenvalues  (Eq.  (2-1))  are  compared  with  the 
experimental  ones  in  the  least-squares  sense.  If  the  differences  in 
eigenvalues 

AX  i = X^(obsd.)  - X^(calcd.) 

are  small,  corrections  of  the  order  of  Af  can  be  applied  to  the  force 
constant  matrix  F to  reduce  the  error  vector  AX  (85) . Statistical 
methods  are  used  to  evaluate  the  results  so  obtained  (26,85-87). 
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Unfortunately,  for  larger  molecules  even  the  frequency  data  are 
incomplete  and  the  additional  data  are  non-existent.  As  a consequence, 
it  is  necessary  to  impose  constraints  on  the  force  field  (12)  which 
are  often  very  difficult  to  justify  physically.  Even  with  these 
additional  constraints,  the  iterative  least-squares  procedure  often 
becomes  unstable  and  the  refinement  fails  to  converge. 

The  principal  problem  with  the  least-squares  method  is  the  fact 
that  it  usually  only  fits  the  fundamental  observed  frequencies  (3N  - 6, 
where  N is  the  number  of  atoms),  which  will  be  subject  to  experimental 
uncertainties.  Moreover,  for  large  molecules  the  additional  experi- 
mental data  (frequencies  from  isotopic  molecules,  Coriolis  constants, 
etc.)  are  not  available  and  it  is  not  possible,  at  this  stage,  to  do 
better . 

Let  us  consider  a second  approach  to  the  problem  of  force  constant 
determination.  This  involves  the  theoretical  evaluation  of  the 
molecular  energy  by  quantum  mechanical  methods  at  the  jab  initio  level 
(76) . Recent  progress  in  computational  techniques  and  computers  has 
made  reliable  ab  initio  force  field  calculations  increasingly  accesible 
to  spectroscopists  (14,40,73).  This  has  significantly  improved  our 
knowledge  of  molecular  force  fields. 

Both  approaches  for  determining  force  constants  will  be  examined 
in  more  detail  in  the  following  discussion.  Some  examples  will  be 
presented,  but  our  principal  concern  will  be  to  examine  the  validity 
of  each  of  these  approaches  for  treating  moderate-to-large  sized 


molecules . 
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Least-square  Refinements  Approach 

In  the  past,  the  usual  method  of  obtaining  the  force  constants 
for  a large  molecule  has  been  to  transfer  values  from  similar  and 
usually  simpler  molecules,  and  then  refine  to  fit  the  observed  fre- 
quencies. A starting  point  would  be,  for  example,  to  transfer  force 
constants  from  small  molecules  with  analogous  fragments,  since  the 
experimental  information  available  for  a simpler  molecule  may  be 
complete  enough  to  determine  accurately  that  molecule's  force  constants. 
These  values  then  form  the  initial  basis  for  a sequence  of  improvements 
or  refinements  generated  by  comparison  with  experimental  data.  There 
is  a vast  literature  related  to  this  approach  (12-14,26,85,89-96).  It 
was  pointed  out  by  Needham  and  Overend  (92)  that  because  the  experi- 
mental information  is  usually  incomplete  for  all  but  the  simplest 
molecules,  the  result  of  a refinement  of  force  constants  contain  a 
large  measure  of  arbitrariness.  In  the  particular  case  where  the 
number  of  independent  pieces  of  data  greatly  exceed  the  number  of 
independent  force  constants  (in  simple  molecular  systems)  this  approach 
has  met  with  considerable  success  (87,93). 

Let  us  consider  some  of  these  cases  where  the  application  of  such 
treatments  has  produced  reliable  force  constants.  For  example,  a refined 
force  field  for  NF^  was  obtained  by  Allan,  Duncan,  Holloway,  and  McKean 
(94)  using  isotope  shifts.  They  were  able  to  rule  out  any  ambig- 
uous force  constants  by  determining  the  centrifugal  distortion  constants 
of  the  molecule.  Their  force  field  reproduced  quite  well  the  experi- 
mental frequencies.  Another  good  example  is  the  study  by  Jones,  Ryan, 

and  Asprey  (95)  to  determine  the  force  constants  of  NOF,  and  NOBr.  They 

18  15 

used  frequency  data,  centrifugal  distortion  constants,  and  0,  N,  and 
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halogen  isotope  shifts  to  fit  the  force  constants  to  the  spectral  data. 

A third  and  final  example  is  the  force  field  determined  for  CC^  by 
Machida  and  Overend  (96) . They  were  successful  in  obtaining  force 
constants  taking  into  consideration  anharmonicity  effects.  A great 
many  other  refined  and  reliable  force  fields  of  small  molecules  (6 
atoms  or  less)  have  been  reported  in  the  literature  for  the  past  forty 
years.  For  more  examples  see  the  reviews  by  Shimanouchi  (86),  Duncan 
(88),  Shimanouchi  and  Nakagawa  (89),  and  Overend  (90). 

Since  the  force  constant  refinement  approach  relies  heavily  upon 
the  observed  frequencies,  the  frequencies  for  a large  number  of  similar 
molecules  (including  isotopic  derivatives)  must  be  measured  and  assigned 
in  order  to  achieve  a meaningful  refinement.  Such  complete  sets  of 
data  are  not  presently  available  for  any  but  a very  few  small  mole- 
cules, and  biological  molecules,  in  particular,  have  not  been  so 
thoroughly  studied. 

A number  of  least-square  force  constants  refinements  for  the 
pyrimidine  (12  atoms)  and  purine  (13  - 16  atoms)  bases  have  been  made 
(1,12,13),  with  resulting  fits  that  are  not  nearly  so  good  as  those  for 
the  simpler  molecules  given  as  examples  above.  For  example,  the 
refinement  of  Susi  and  Ard  (12)  for  uracil  gives  a rather  poor  fit  in 
the  frequency  region  900  cm  ^ to  1400  cm  ^ . Furthermore,  they  found  it 
necessary  to  constrain  many  of  the  force  constants  to  constant  values 
during  the  refinement  in  order  to  prevent  the  refinement  from 
diverging  as  the  iteration  proceeded. 

Quantum  Mechanical  Approach 

Force  constants  are  determined  by  the  electronic  nature  of  the 
atoms  forming  the  bonds  in  a molecule.  There  is  some  reason  to  believe 
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that  chemically  similar  parts  of  molecules  (for  example  a C-H  bond  in 
ethane  and  a C-H  bond  in  methane)  should  have  similar  force  constants. 

The  ab  initio  quantum  mechanical  calculation  of  force  constants  in 
recent  years  has  supported  this  idea  (40,73,75,97).  For  the  present 
work,  it  is  of  special  significance  that  a force  field  obtained  for  a 
molecule  in  an  ab  initio  calculation  using  the  4-3 1G  basis  set  reproduces 
quite  well  the  spectral  frequencies,  when  multiplied  by  the  proper 
scaling  factors  (14,98-100).  Pulay,  Fogarasi,  Pongor,  and  Vargha  (97) 
have  recommended  a uniform  procedure  to  carry  out  such  scaling  proce- 
dure for  results  from  a similar  4-21G  basis  set.  Their  determination 
of  the  scaling  factors  is  based  on  a least-squares  fitting  of  the 
calculated  vibrational  frequencies  to  the  experimental  ones.  Empiri- 
cally adjusted  ab  initio  force  constants  are  now  available  for  a number 
of  prototype  organic  molecules  such  as  benzene  (73),  methanol  (98), 
propene  (99),  ethane  (100),  cyclobutane  (101),  and  cubane  (102).  It  is 
concluded  that  a combination  of  theoretical  calculations  with  experi- 
mentally assigned  frequencies  makes  it  possible  to  obtain  significant 
values  of  adjusted  force  constants.  For  example,  it  is  generally  found 
that  diagonal  stretching  force  constants  are  overestimated  in  ab  initio 
calculations  by  about  10  - 15  % and  diagonal  bending  force  constants 
are  overestimated  by  about  20  - 30  %,  but  occassionally  calculated 
interaction  force  constants  agree  very  well  with  the  experimental 
values.  These  adjusted  ab  initio  force  fields  have  proved  to  be  very 
useful  for  predicting  spectra  of  molecules  (73,74,98). 

It  should  be  made  clear  that  the  scaling  procedure  for  the 
ab  initio  results  can  be  carried  out  either  for  force  constants  or  for 


the  frequencies.  For  various  reasons,  the  former  approach  is  commonly 
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preferred.  First,  frequencies  can  be  determined  from  the  force 
constants  but  not  vice  versa.  Also,  many  different  types  of  force 
constants  simultaneously  contribute  to  a specific  frequency.  Moreover, 
the  "correction"  of  ab  initio  force  constants  with  the  accuracy  given 
for  the  scaling  factors  in  the  preceding  paragraph  still  leads  to  over- 
estimation of  the  vibrational  frequencies  by  roughly  10  %.  Pople, 
Schlegel,  Krishnan,  De  Frees,  and  Binkley  (103)  demonstrated  this 
conclusion  in  a study  of  the  calculated  frequencies  of  37  molecules. 
Other  studies  (40,74,104)  have  reached  the  same  conclusion. 

So  far  we  have  seen  that  the  scaling  factors  must  play  a major 
role  in  any  scheme  involving  the  transfer  of  "adjusted  force  constants" 
of  smaller  molecules  to  predict  vibrational  frequencies  of  larger  mole- 
cules. To  illustrate  further  we  present  the  following  sample  of 
studies  which  have  adjusted  force  constants  obtained  from  ab  initio 
quantum  mechanical  methods  by  empirical  scale  factors. 

Scaled  ab  initio  force  fields  were  systematically  obtained  by  Blom 
and  Altona  for  ethane,  propane,  cyclopropane  (100),  propene  (99), 
methanol  (98),  and  ethene  ( 105) . In  all  their  calculations  they  used 
the  4-31G  basis  set  to  obtain  the  complete  general  valence  force  field 
of  each  molecule.  This  ab  initio  force  field  was  later  scaled  with 
different  scaling  factors  in  order  to  reproduce  the  experimental 
frequencies.  For  the  first  three  molecules  (100)  a total  of  180 
vibrational  frequencies  (including  those  for  deuterated  analogues) 
were  used  in  a fitting;  the  mean  deviation  between  observed  and  cal- 
culated frequencies  amounted  to  10.4  cm  1 or  0.74  %.  Similar  results 


were  obtained  for  the  other  molecules. 
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For  the  benzene  molecule,  Pulay,  Fogarasi,  and  Boggs  (73)  obtained 
ab  initio  force  fields  with  different  basis  sets  (ST0-3G,  4-21G,  and 
4-31G) , which  were  then  scaled  with  scaling  factors  determined  by 
least-squares  fitting  to  the  observed  spectrum.  They  found  that  these 
empirical  scale  factors  fall  in  the  range  0.7  — 0.9.  Another  important 
conclusion  from  this  study  was  the  fact  that  the  calculated  ST0-3G 
force  constants  are  qualitatively  correct,  in  that  at  least  their  signs 
agree  with  those  from  the  better  basis  sets. 

The  scale  factors  obtained  for  benzene  (73)  were  directly  applied 
to  the  ab  initio  (4-21G  basis  set)  force  field  constants  for  pyridine 
by  Pongor,  Pulay,  Fogarasi,  and  Boggs  (106).  This  scaled  force  field 
was  then  used  to  predict  the  vibrational  spectrum  of  pyridine.  Compar- 
ison of  experimental  and  predicted  spectra  showed  a mean  deviation  of 
5.7  cm  for  the  non-CH  frequencies.  They  suggested  that  this  method 
of  obtaining  scale  factors  for  a related  reference  molecule  and 
applying  them  to  another  will  permit  the  accurate  prediction  of  a truly 
unknown  vibrational  spectrum. 

Using  an  .ah.  initio  calculation  with  a 4-31G  basis  set,  Wiberg, 
Walters,  Wong,  and  Colson  (74)  also  determined  a force  field  for 
pyridine,  which  they  adjusted  with  scaling  factors  (constant  for  each 
symmetry  class,  and  ranging  from  0.73  - 0.86)  to  fit  the  observed 
frequencies  for  pyridine-d^  , pyridine-d^,  and  partially  deuterated 
pyridines.  Their  fit  between  the  observed  and  calculated  frequencies 
gave  a standard  deviation  of  6.6  cm-'*'. 

Transfer  of  Atomic  Polar  Tensors 

As  it  has  been  previously  pointed  out  in  Chapters  2 and  3,  the  APTs 
for  the  atoms  in  a molecule  are  a natural  representation  of  intensity 
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results  for  comparison  with  ab  initio  quantum  mechanical  calculations 
for  that  particular  molecule.  Since  we  are  concerned  here  with  the 
possibility  of  successful  transfer  of  vibrational  parameters,  we  shall 
address  the  following  question:  should  we  expect  that  the  APTs  for  the 
atoms  in  a molecule  to  be  exactly  transferable  to  the  atoms  in  a 
chemically  similar  molecule?  If  the  evidence  demonstrates  that  the 
answer  to  this  question  is  affirmative,  then  the  transferability  of 
APTs  provides  a powerful  tool  for  the  prediction  of  absolute  infrared 
intensities  of  larger  molecules,  starting  with  APT  values  from  smaller 
molecules.  As  we  attempt  to  answer  this  question  it  should  be  remem- 
bered that  the  elements  of  the  APTs  are  expected  to  reflect  contribu- 
tions from  dipole-induced  dipole  terms,  changes  in  hybridization, 
etc.,  which  are  not  all  expected  necessarily  to  be  the  same  from  one 
molecule  to  another. 

Transfer  of  APTs  from  reference  molecules  with  the  aim  of 
calculating  the  gas  phase  infrared  intensities  for  similar  molecules 
has  been  discussed  by  Person  (60) . He  pointed  out  that  since  general 
rules  have  not  been  established  to  decide  which  APTs  are  the  most 
suitable  to  be  transferred  to  the  molecule  whose  intensities  are  to  be 
determined,  this  question  poses  a real  challenge.  However,  some 
attempts  by  person  and  coworkers  (23,32,39,63,64,107),  Bruns  (108,109), 
and  others  (58)  have  provided  promissing  results  which  might  well 
serve  as  guidelines  when  trying  to  transfer  APTs  of  atoms  in  a molecule 
to  the  same  atoms  in  another  molecule. 

The  first  suggestion  that  APTs  may  be  transferable  was  the  dis- 
covery by  King,  Mast,  and  Blanchette  (62)  that  the  effective  charge, 

£ , appeared  to  be  the  same  for  H atoms  in  a wide  variety  of  molecules. 
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Since  the  E,a  (later  redefined  in  terms  of  xa  as  three  times  xa>  (see 
(3-6)  and  discussion  of  APT  invariants  in  Chapter  3))  is  an  invariant  of 
the  APT,  Newton  and  Person  (32)  conducted  a study  to  examine  the  possi- 
bility of  transferring  the  APTs  themselves  for  H and  F atoms  in  differ- 
ent fluorinated  methane  molecules.  Their  results  were  very  encouraging. 
They  concluded  that  the  APTs  for  H and  for  F atoms  from  CtT^F  are  reason- 
able choices  for  transfer  to  calculate  the  infrared  intensities  of  other 
f luoromethane  molecules.  They  calculated  intensities  for  CH^F^j  CHF^, 
and  CF^  using  the  APTs  from  CI^F.  The  results  obtained  were  always 
within  a "factor  of  two"  (32,63)  of  the  experimental  values,  and  for 
strong  bands  the  agreement  is  much  better. 

Additional  studies  by  Person  and  Newton  (64,110)  showed  that  this 
F atom  APT  could  be  transferred  to  predict  the  infrared  intensities  for 
fundamental  vibrations  of  the  CF2  and  CF^  radicals.  It  was  not  obvious 
that  the  use  of  the  F atom  APT  in  this  calculation  would  be  successful, 
since  CF^  has  an  odd  electron,  and  CF^  has  two  odd  electrons.  Never- 
theless, the  results  obtained  are  in  very  good  agreement  with  exper- 
imental data.  It  has  been  pointed  out  by  Person  (60)  that  this  study 
illustrates  the  potential  power  of  this  technique  for  applications  to 
infrared  studies  of  free  radical  reactions.  Furthermore,  the  F atom 
APT  of  CH^F  was  transferred  to  predict  the  infrared  intensities  for 

SF,  and  UF,  molecules,  by  Person  and  Overend  (111).  They  found  that 
6 6 

their  predicted  intensities  agreed  with  experimental  values  within  the 
experimental  uncertainty. 

Now  let  us  consider  some  examples  where  the  F atom  APT  from  CH  F 
cannot  be  successfully  transferred  to  make  intensity  predictions.  This 
is  the  case  for  the  CF^S,  molecule,  as  has  been  demonstrated  by  Golden, 
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Horner,  and  Overend  (112).  This  non-transferability  has  been  discussed 
by  Person  (60)  in  terms  of  atomic  polarizabilities.  He  points  out  that 
we  should  not  expect  success  in  the  transfer  of  an  F atom  APT  from  CH^F 
to  CF3C£  because  C t is  more  polarizable  than  H or  F atoms,  so  this 
"polarizability  effect"  causes  a big  charge  flux  contribution  to  the 
F atom  APTs  in  CF3C£  but  not  in  C^F  or  in  CF^  The  same  argument  has 
been  used  to  explain  failures  of  transferability  of  APTs  to  other 
molecules;  for  example,  in  F^CO  and  F^CS  (113). 

It  has  been  shown  by  Rogers  (39)  that  H and  0 atom  APTs  are  also 

transferable  from  CH  F and  ^0  to  C^OH,  C^OH,  and  CH3OCH3  molecules. 

The  transfer  was  not  successful  from  HCHO  to  CH3COCH3,  probably  because 

of  the  influence  from  the  polarizable  C=0  bond  on  the  C-H  intensity 

parameters.  It  was  also  found  by  Herrera,  Ramos,  and  Bruns  (108)  that 

the  APTs  from  cis-dif luoro—  and  cis-dichloroethylene  are  directly 

transferable  to  their  respective  trans  isomers,  resulting  in  predicted 

intensities  that  are  in  agreement  with  the  experimental  intensities. 

In  another  study,  Kagel,  Powell,  Hopper,  Overend,  Ramos,  Bassi,  and 

Bruns  (114)  successfully  transferred  H,  F,  and  0 atom  APTs  to  HFCO  and 

DFCO  molecules  from  cis-C^H^F^,  ^CO  and  F^CO  molecules.  They  reported 

that  the  calculated  values  of  the  infrared  intensities  for  the  former 

molecules  are  in  good  agreement  with  the  experimental  values.  Also, 

transferring  the  C and  0 atom  APTs  of  Cr(C0)6  for  the  prediction  of  the 

intensities  of  Mo (C0)£  and  W(C0),,  Bruns,  Hase,  and  Brin  (109)  repro— 

o o 

duced  the  experimental  intensities  of  these  heavier  group  VIB  metal 
carbonyls  within  the  usual  "factor  of  two".  However,  when  they 
attempted  a direct  transfer  of  the  APTs  from  Ni(CO)^  to  Mo(CO)^  and 

W(C0) , the  experimental  intensities  were  not  correctly  reproduced. 
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They  suggested  first,  that  this  failure  could  be  attributed  to  the 
difference  in  electronic  structure  of  the  Ni  carbonyl  as  compared  to 
that  of  Cr,  Mo,  or  W carbonyls;  and  second,  that  the  different 
symmetries  imply  differences  in  hybridization  of  the  metal  orbitals, 
and  therefore  different  behaviours  with  respect  to  variation  of  bond 
distances  and  angles. 

All  the  evidence  presented  shows  that  a.  priori  guidelines  have  not 
been  established  as  yet  for  transferring  the  appropriate  APTs  for 
"perfect"  intensity  predictions.  However,  it  is  very  encouraging  to 
note  that  the  evidence  demonstrates  succesful  results  within  a "factor 
of  two"  in  the  prediction  of  infrared  intensities  when  the  APTs  are 
transferred  for  atoms  in  chemically  related  molecular  environments. 

Comparison  of  ab  initio  Vibrational  Parameters 

An  examination  and  comparison  of  the  force  constants  and  APTs 
obtained  by  ab  initio  quantum  mechanical  methods  for  uracil  and  other 
molecules  (e.g.,  pyridine)  is  expected  to  provide  interesting  details 
about  the  similarity  or  difference  in  the  electronic  nature  of  these 
molecules.  Perhaps  we  could  also  obtain  information  on  the  direct 
transfer  of  such  vibrational  parameters  between  related  molecules.  Of 
course,  we  are  not  looking  for  exact  values  of  ab  initio  force 
constants  or  APTs  in  two  molecules,  but  rather  for  possible  trends  in 
values  of  these  parameters.  Anyway,  this  type  of  comparison  should 
give  more  confidence  about  the  use  of  quantum  mechanical  procedures  in 
solving  vibrational  spectroscopy  problems. 

Ab  initio  Force  Constants 

The  first  comparison  of  ab  initio  values  that  we  want  to  make  is 
for  some  "corrected  to  4— 31G"  force  constants  (see  Chapter  4)  of  uracil 
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with  those  obtained  for  other  molecules  at  the  same  level  of  approxi- 
mation. In  Table  5-1  are  presented  and  compared  some  of  these  force 
constants.  The  values  that  are  given  in  each  column  were  taken  from 
the  following  studies,  Column  1:  ST0-3G  "converted  to  4-31G"  for  uracil 
by  Nishimura,  Tsuboi,  Kato,  and  Morokuma  (14);  Column  2:  4-31G  force 
constants  determined  for  formamide  by  Sugawara,  Hamada,  Hirakawa, 
Tsuboi,  Kato,  and  Morokuma  (46);  Column  3:  4-31G  values  for  N-methyl- 
formamide  obtained  by  Sugawara,  Hirakawa,  Tsuboi,  Kato,  and  Morokuma 
( 4 0) ; Column  4:  gives  some  force  constants  at  the  4-31G  level  for 
pyridine  by  Wiberg,  Walters,  Wong,  and  Colson  (74);  Column  5:  shows 
some  values  at  the  4-31G  level  for  acetaldehyde  obtained  by  Wiberg, 
Walters,  and  Colson  (104);  and  Column  6:  gives  some  force  constants 
obtained  for  pyridine  at  the  4-21G  level  by  Pulay,  Fogarasi,  and  Boggs 
(106). 

A close  look  to  Table  5-1  shows  how  similar  the  value  of  the  ab 
initio  force  constants  are  in  the  different  molecules  presented.  Take, 
for  example,  the  C^O  stretching  values  of  13.298  md  A ^ for  uracil, 
13.864  md  A ^ for  formamide,  13.4  md  A * for  N-methylformamide,  and 
14.749  md  A * for  acetaldehyde.  It  is  evident  that  all  these  values 
for  this  constant  fall  in  a very  narrow  range.  The  same  similarity 
holds  true  for  the  force  constants  for  the  C2N^  stretch,  N^H  defor- 
mation, and  C^N  deformation.  It  is  also  interesting  to  note  that  the 
"corrected  to  4-31G"  force  constant  value  of  6-.  181  md  A ^ for  the  C^H 
stretch  in  uracil  seems  to  be  overestimated  when  compared  to  the  rest 
of  the  values  given  in  Table  5-1:  5.708,  5.77,  5.218,  5.198,  and 
5.658  md  A ^ . 


TABLE  5-1 

COMPARISON  OF  AB  INITIO  FORCE  CONSTANTS 
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In  general,  the  interaction  force  constants  calculated  for  the 
various  molecules  presented  in  Table  5-1  also  have  comparable  values 
within  the  same  order  of  magnitude.  Also,  their  signs  are  usually- 
predicted  the  same  in  all  cases.  It  is  not  clear  whether  or  not  the 
discrepancies  are  meaningful. 

APTs  for  uracil  and  pyridine 

Let  us  now  compare  the  APTs  for  uracil  (see  Table  3-3)  from  our 
calculations  (17)  with  those  for  pyridine  obtained  by  Wiberg,  Walters, 
Wong,  and  Colson  (74).  Both  calculations  were  carried  out  using  a 

4- 3 1G  basis  set.  In  table  5-2  the  APTs  for  each  atom  in  uracil  and 
pyridine  are  shown.  The  atom  numbering  follows  that  given  in  Fig. 

5- l(a).  In  order  to  facilitate  their  comparison,  the  APTs  have  been 
expressed  in  a bond  coordinate  system,  as  defined  in  Figs.  5-l(b-g) . 
Also,  the  APTs  for  pyridine  obtained  by  Wiberg  et  al.  (74)  have  been 
converted  from  D £ ^ to  e units  (le  = 4.803  D X S- 

We  see  in  Table  5-2  that  the  values  for  the  APTs  are  very  differ- 
ent for  almost  all  atoms.  It  seems  that  the  aromaticity  of  pyridine 
and  non-aromaticity  of  uracil  causes  the  atoms  in  the  ring  (N^ , C2>  C^, 
C^,  and  C^)  to  be  different  chemically  in  the  two  rings  and  thus  to 
have  different  effective  charges  x-  Basically  the  only  similarities 
between  these  molecules  lie  in  the  APTs  and  the  x values  of  atoms 
and  H^2*  This  could  well  be  expected  given  that  the  double  bond 

in  uracil  could  resemble  the  analogous  bond  in  pyridine  in  its  elec- 
tronic nature  and  polarizability  character.  Hence,  the  terminal 
hydrogen  atoms  attached  to  these  carbon  atoms  should  be  very  similar. 

It  is  interesting  to  notice  that  the  atom  is  so  different  in 
pyridine  than  in  uracil,  due  partly  to  attached  H-,  atom  in  the  latter. 


TABLE  5-2 

COMPARISON  OF  APTs  FOR  URACIL 
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Figure  5-1.  Coordinate  bond  axes  representation,  (a)  Atom 
numbering  for  uracil  and  pyridine,  (b)  Bond 
coordinate  system  for  atoms  2 and  8.  Also, 
(c),  (d),  (e),  (f),  and  (g)  represent 
bond  coordinate  systems. 
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As  shown  in  Table  5-2,  a sum  of  the  APT  elements  of  atoms  ^ and  H7  of 
uracil  gives  some  values  similar  to  those  of  the  APT  for  pyridine.  In 
particular  the  ( 3p  /9x)  and  (9p  /9z)  values  of  -0.488  and  -0.388  e, 
respectively,  for  pyridine  compare  very  well  with  the  sum  of  the 
corresponding  elements  for  atoms  and  in  uracil  (—0.480  and  —0.292 
e,  respectively) . Also  note  that  the  effective  charge  x for  atom  ^ in 
uracil  is  2.16  times  that  of  the  N atom  of  pyridine,  while  the  ratio 
of  x(Nx  + h7  for  uracil)  to  x(N1  for  pyridine)  is  just  1.66. 

It  is  also  evident  from  Table  5—2  that  the  APTs  of  the  carbon 
atoms  on  the  ring  of  both  molecules  show  many  fluctuations.  The 
values  for  these  atoms  in  uracil  are  always  larger  than  those  in 
pyridine. 

In  conclusion,  from  this  comparison  is  obvious  that  the  transfer 
of  pyridine  parameters  to  a pyrimidine  ring  would  produce  incorrect 
results. 

Transfer  of  ab  initio  Vibrational  Parameters 
to  Methylated  Uracifs 

Because  of  our  interest  in  the  methylated  uracil  molecules,  we 
shall  end  this  chapter  discussing  how  force  constants  and  APTs  might  be 
transferred  to  predict  spectra  of  these  particular  molecules.  This 
survey  will  serve  as  the  basis  for  presentation  of  the  actual  results 
in  Chapter  6. 

As  described  in  Chapter  2 we  decided  to  carry  out  two  types  of 
calculations,  namely:  Calculation  I where  the  methyl  group  is  taken 
as  a dynamic  unit  of  mass  15  u,  and  Calculation  II  where  the  methyl 
group  is  properly  treated  as  a four  atom  dynamic  unit. 
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Transfer  of  ab  initio  Force  Constants 

The  methylated  uracils  to  be  considered  are  the  following: 
1-methyluracil  ( 1-Me  U) , l-methyl-3-deuterouracil  (l-Me-3-D  U) , 

3-methyl uracil  (3-Me  U)  , and  1 , 3-dimethyluracil  (1,3-di  Me  U) . Each 
one  of  these  molecules  will  be  discussed  in  more  detail  in  Chapter  6. 
Here  we  just  want  to  discuss  which  force  constants  can  be  used  in 
each  case. 

For  1-Me  U we  can  transfer  the  uracil  "corrected  to  4-31G"  force 
constants  directly  to  the  ring  system  of  the  molecule,  when  carrying 
out  both  Calculation  I and  Calculation  II.  The  same  transfer  should 
apply  to  the  rings  of  1-Me- 3-D  U,  3-Me  U,  and  1,3-di  Me  U.  When 
treating  any  of  these  molecules  in  Calculation  I,  we  can  transfer  the 
stretching,  bending,  and  interaction  ab  initio  force  constants  from 
N-methylformamide  (40,  see  also  Table  5-1)  to  the  corresponding  N-Me 
bond  in  the  methylated  uracils. 

The  values  of  the  C-H  stretching,  bending,  and  interaction  force 
constants  to  be  used  for  the  methyl  group  in  Calculation  II  can  be 
directly  transferred  from  the  CH^  group  attached  to  the  nitrogen  atom 
in  N-methylformamide,  from  the  ab  initio  calculation  of  the  force 
constant  values  for  this  molecule  (40) , obtained  using  a 4-31G  basis 
set.  Also,  the  comparison  of  ab  initio  force  constants  in  Table  5-1 
shows  some  consistency  between  the  values  of  force  constants  for  uracil 
and  for  N-methylformamide.  Another  possibility  is  to  make  a direct 
transfer  of  the  force  constants  for  the  methyl  group  obtained  at  the 
4-31G  level  for  N-methylacetamide  by  Sugawara,  Hirakawa,  Tsuboi,  Kato, 
and  Morokuma  (115) . We  decided  to  use  the  force  constant  values  from 
the  N-methylacetamide  study  (115)  because  we  had  obtained  the  APTs  for 
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N-methylacetamide,  and  we  wanted  to  conduct  a full  study  of  frequencies 
and  intensities  for  this  molecule,  which  will  be  presented  elsewhere 
(116). 

Transfer  of  APTs 

In  order  to  decide  on  which  APTs  should  be  transferred  to  the 
atoms  in  1-Me  U,  l-Me-3-D  U,  3-Me  U,  and  1,3-di  Me  U,  the  simplest 
approach  would  be  just  the  direct  transfer  of  APTs  from  each  atom  in 
uracil  to  the  corresponding  atoms  in  each  particular  molecule.  For 
Calculation  I this  can  be  done  exactly.  The  APT  of  the  NjH  hydrogen 
atom  in  uracil  is  transferred  to  the  Me  dynamic  unit  in  1-Me  U,  or  that 
for  N^H  to  3-Me  U,  etc.  This  assumption  is  perhaps  the  most  drastic  in 
Calculation  I. 

For  Calculation  II  there  is  a bigger  challenge  in  deciding  on  the 

correct  transfer  of  APTs  for  the  N-CH^  fragment  in  each  one  of  the 

methylated  uracils.  The  APTs  for  the  CH^  groups  in  acetone  obtained 

by  Rogers  (39)  from  an  ab  initio  calculation  using  a 4-31G  basis  set 

seem  to  reproduce  well  the  intensity  pattern  of  the  methyl  groups  in 

the  C-H  stretching  region  of  the  methylated  uracils.  We  can  satisfy 

the  null  condition  (Eq.  (3-5))  by  adjusting  the  APTs  of  the  carbon  and 

nitrogen  atoms.  In  other  words,  we  can  transfer  the  APTs  for  all  four 

atoms  in  the  complete  CH^  group  from  acetone  to  the  methylated  uracils 

and  use  the  null  condition  to  determine  the  APT  of  the  nitrogen  atom  to 

which  the  CH  is  attached,  or  we  might  transfer  those  APTs  for  the 
3 

hydrogen  atoms  from  acetone  to  the  methylated  uracils,  using  the  null 
condition  to  obtain  the  APT  for  the  carbon  atom  of  the  CH^  group.  We 
have  considered  both  alternatives  in  this  study.  Again  we  emphasize 
that  the  APTs  for  the  rest  of  the  atoms  are  transferred  directly  from 


uracil . 
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Inspired  by  the  ab  initio  study  of  force  constants  for  N-methyl- 
formamide  (40)  we  decided  to  conduct  an  ab  initio  study  of  the  APTs  of 
a related  molecule,  N-methylacetamide  (116).  Since  this  molecule 
(N-methylacetamide)  has  two  CH^  groups , we  propose  to  test  the  transfer 
of  APTs  from  each  one  of  them  to  the  CH^  group  in  the  methylated 
uracils.  Again,  the  null  condition  can  be  used  to  determine  the 
nitrogen  and  carbon  APTs  in  the  N-CH^  moieties  of  the  methylated 
uracils.  The  APTs  for  N-methylacetamide  are  given  in  Table  5-3,  from 
Ref.  116.  The  atom  numbering  and  axes  orientation  for  the  APTs  in 
Table  5-3  is  according  to  Fig.  5-2. 


TABLE  5-3 

APTs  FOR  N-METHYLACETAMIDE 
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Note:  The  atomic  numbering  scheme  and  coordinate  system  are  shown  in  Fig.  5-2.  Effective  charges  (x)  are 
calculated  from  Eq.  (3-6).  Units  are  electrons  (e). 
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Figure  5-2.  Atom  numbering  and  orientation  of  the  principal 
Cartesian  coordinate  axes  of  N-methylacetamide . 


CHAPTER  6 

PREDICTION  OF  INFRARED  SPECTRA  OF  METHYLATED  URACILS 

This  chapter  deals  with  the  details  and  results  of  the  prediction 
of  infrared  spectra  for  the  methylated  uracils,  using  both  Calculations 
I and  II.  We  shall  first  consider  the  transfer  of  APIs  and  ab  initio 
force  constants  from  smaller  molecules  or  molecular  fragments  to  the 
following  methylated  uracils:  1-methyluracil  (1-Me  U) , l-methyl-3- 
deuterouracil  (l-Me-3-D  U)  , 3-methyluracil  (3-Me  U) , and  1,3-dimeth- 
yluracil  (1,3-di  Me  U) . This  is  followed  by  a discussion  and  compar- 
ison of  the  results  of  the  experimental  matrix-isolated  and  predicted 
spectra  for  these  molecules,  in  order  to  assign  their  observed  frequen- 
cies and  intensities. 

Transfer  of  APTs 

In  our  scheme  for  transferring  vibrational  parameters  presented  in 
Chapter  2 we  explained  that  our  predictions  of  frequencies  and  intensi- 
ties were  made  on  two  levels:  Calculation  I and  Calculation  II.  Let  us 
now  consider  the  transfer  of  APts  to  the  methylated  uracils  using  each 
one  of  these  approaches. 

Calculation  I 

In  the  first  and  simpler  approach  (Calculation  I)  the  methyl  group 
was  treated  as  a single  dynamical  point  mass  of  15  u,  so  this  calcula- 
tion gives  the  effect  of  the  heavy  mass  on  the  frequencies  of  the  ring 
vibrations.  Basically,  the  problem  to  solve  is  one  similar  to  that  of 
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determining  frequencies  and  intensities  for  a deuterated  uracil  molecule 
except  that  we  use  15  u for  the  substituted  atom  instead  of  2 u.  In  the 
experimental  spectrum  of  the  real  molecule,  of  course,  we  see  the  charac 
teristic  absorptions  by  the  vibrations  within  the  methyl  group,  but  we 
assume  that  this  spectrum  is  known  because  it  is  characteristic.  After 
making  the  changes  in  masses  in  uracil  (Fig.  3-1),  we  place  an  atom  of 
15  u in  the  Hy  position  for  1-Me  U and  l-Me-3-D  U (in  this  molecule  2 u 
is  the  mass  for  the  replacement  of  Hg  by  Dg) , in  the  Hg  position  for 
3-Me  U,  and  two  15  u atoms  replace  Hy  and  Hg  for  1,3-di  Me  U.  Basical- 
ly, the  force  constants  are  assumed  to  be  the  same  as  uracil,  and  the 
frequencies  are  calculated  for  the  "heavy  mass  uracils  . 

All  the  APTs  used  in  Calculation  I to  predict  intensities  were 
transferred  directly  from  uracil  to  the  methylated  uracils.  These  APTs 
are  shown  in  Table  3-3  for  uracil.  No  rotation  of  these  APTs  was  ne- 
cessary, because  we  made  our  calculation  in  the  same  Cartesian  coordi- 
nate system  which  was  used  for  uracil  (Fig.  3-1). 

Calculation  II 

The  approach  followed  in  Calculation  II  is  much  more  elaborated 
because  the  methyl  group  was  treated  properly  as  a polyatomic  substi- 
tuent, so  we  have  to  carry  out  NCAs  for  15-atom  or  18-atom  molecules 
instead  of  for  12-atom  molecules.  The  APTs  used  for  the  CH3  group  in 
this  calculation  were  transferred  from  acetone  (39)  or  N-methyl- 
acetamide  (116),  and  those  for  the  ring  atoms  were  transferred  from 
uracil . 

1-methyluracil 

Let  us  consider  first  the  case  of  Calculation  II  for  1-Me  U,  in 
which  Hy  in  uracil  (Fig.  3-1)  is  replaced  by  the  CH3  group.  The 
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equilibrium  geometry  and  the  numbering  of  atoms  for  this  molecule  are 
shown  in  Fig.  6-1,  where  the  structural  parameters  given  for  the  N-CH3 
fragment  are  derived  from  those  reported  for  that  fragment  in  N-methyl- 
formamide  by  Kitano  and  Kuchitsu  (117)  from  an  electron  diffraction 
study.  Khetrapal  and  Kunwar  (118)  found  from  NMR  measurements  that  the 
I-CH3  group  in  1,3-di  Me  U is  oriented  with  one  hydrogen  in  the  plane 
pointing  toward  the  2-carbonyl  group,  and  we  have  assumed  this  geometry 
for  1-Me  U.  Also  included  in  Fig.  6-1  are  the  Cartesian  coordinates  for 
the  CH3  group;  the  coordinates  for  the  remaining  atoms  are  the  same  as 
for  uracil  (Fig.  3-1  and  Table  3-1) . 

In  order  to  test  the  transferability  of  APTs  for  the  methyl  in  the 
methylated  uracils  in  Calculation  II,  we  tried  several  different  sets  of 
APTs,  changing  particularly  those  for  the  atoms  in  the  CH3  group  from 
the  values  transferred  from  acetone,  to  those  from  the  N-CH3  and  from 
the  C-CH3  groups  of  N-methylacetamide  (Table  5-3) . For  1-Me  U the  APTs 
transferred  from  acetone  were  those  calculated  by  Rogers  (39) , but  the 
coordinate  system  was  rotated  by  +87.74°  (Eq.  (3-4))  in  order  to  have 
them  coincide  with  the  coordinate  axes  for  1-Me  U,  as  seen  in  Fig.  6-2. 
Since  the  APTs  for  N-methylacetamide  (116)  were  calculated  in  the  prin- 
cipal Cartesian  system  of  Fig.  5-2  (Table  5-3),  it  was  necessary  to 
rotate  this  coordinate  system  for  the  APTs  in  the  N-CH3  group  by 
+232.91°  and  for  the  APTs  in  the  C-CH3  group  by  +60.21°  (Eq.  (3-4)), 
as  it  is  also  shown  in  Fig.  6-2,  in  order  to  conform  the  coordinate 
axes  orientation  of  1-Me  U (see  Fig.  6-1).  The  resulting  APTs  after 
the  rotations  are  given  in  Table  6-1,  where  Column  1 shows  the  APTs 
from  acetone,  and  Columns  2 and  3 those  from  the  N-CH3  and  C-CH3 
groups  of  N-methylacetamide,  respectively.  Also,  the  APT  for  the  N3 
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Figure  6-1.  Geometrical  parameters  of  1-Me  U and  Cartesian  coordinates 
for  the  CH3  group  in  the  Cartesian  coordinate  system  of 
uracil.  Coordinates  for  the  ring  atoms  of  uracil  are  given 
in  Table  3-1.  All  C-H  bond  lenghts  in  the  CH3  group  are 
taken  as  1.092  ft. 
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TABLE  6-1 

TRANSFERRED  APTs  (in  e)  FROM  THE  CHo  GROUPS  OF  ACETONE, 
AND  FROM  THE  N-CH3  AND  C-CH3  GROUPS  OF  N-METHYLACETAMIDE , 
AND  THAT  FOR  THE  N.  ATOM  OF^URACIL,  USED  FOR  THE  N-CH3 
GROUP  OF  1-METHYLURACIL 
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0 ^ 
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0 \ 

-.055 

-.049 

0 

p*<Hn| 

-.102 

-.044 

0 

-.092 

0.042 

0 

-.153 

-.052 

0 

1 

\ 0 

0 

0.101, 

l o 

0 

0.126, 

\ 0 

0 

0.073 

X 

= 0.094 

X 

= 0.101 

X 

= 0.111 

J0.O13 

0.004 

o \ 

-.042 

0.080 

o \ 

<0.057 

-.003 

0 

p*<Hu: 

0.016 

0.081 

0.013 

0.080 

0.034 

0.017 

0.033 

0.064 

0.019 

y.086 

-.038 

-.028J 

y-  .086 

-.018 

-.080) 

p .061 

-.013 

0.012 

X 

= 0.102 

X 

= 0.128 

X 

= 0.083 

/ Nl* 

uracil 

Nls  N 

-Me  A 

. Cl> 

N-Me  A 

-.957 

0.163 

o \ 

-1.426 

0.461 

0 \ 

vO 

C5> 

o 

0.092 

0 

PK 

(Nj.Cj) 

-.678 

-.978 

0 

0.146 

-.512 

0 

0.063 

1.033 

0 

\ o 

0 

-.729 

\ o 

0 

-.085 

\ ° 

0 

0.5Uj 

X 

= 0.982 

X 

= 1.032 

X 

= 1.314 

aAll  APTs  have  been  been  rotated  to  coordinate  axes  of  Fig.  6-1. 

bAPTs  taken  from  the  study  of  acetone  by  Rogers  (39) . 

CAPTs  taken  from  the  study  of  N-methylacetamide  by  Scott  (116). 

dPx(Hi5)  is  obtained  by  reversing  the  signs  of  the  off-diagonal 
components  of  ^(H^)  involving  z. 
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atom  of  uracil  (Fig.  3-1  and  Table  3-3)  is  included  in  Table  6-1,  and 
compared  with  those  for  the  atom  in  the  N-CH^ , and  to  the  atom  in 
the  C-CH3  fragments  of  N-methylacetamide . Note  that  these  APTs  depend 
also  on  the  nature  of  the  other  groups  attached  to  the  (or  C^)  atom, 
as  well  as  in  the  CH3  group.  For  all  the  APTs  the  effective  charges 
(Eq.  (3-6))  for  each  atom  are  also  given  in  this  table. 

In  our  Calculation  II  we  assumed  that  l-Me-3-D  U has  the  same 
geometry  and  vibrational  parameters  (force  constants  and  APTs)  as 
1-Me  U.  Hence,  the  calculation  of  both  molecules  uses  the  same  APTs 
and  force  constants. 

Let  us  now  compare  the  different  sets  of  APTs  given  in  Table  6-1. 
First,  we  note  that  there  is  some  similarity  in  the  diagonal  elements 
of  the  APTs  for  the  N^  atoms  in  uracil  and  in  N-methylacetamide: 

-0.957,  -0.978,  and  -0.729  e in  uracil  compared  with  -1.426,  -0.512, 
and  -0.815  e respectively  in  N-methylacetamide,  with  x values  of  0.982 
and  1.032  e,  respectively.  When  we  compare  the  APTs  for  in  all  three 
cases  in  Table  6-1,  we  see  how  different  the  x values  are:  0.183  e (ace- 
tone), 0.437  e (N-CH3),  and  0.212  e (C-CH3) . The  elements  of  the  APTs 
for  this  atom  are  different  for  each  case.  However,  note  that  those  for 
acetone,  where  the  CH3  group  is  attached  to  a -C0CH3  fragment,  are  more 
similar  to  those  for  C-CH3  in  N-methylacetamide  (where  the  CH3  is 
attached  to  the  -C0NHCH3  fragment).  The  x values  are  quite  close,  in 
fact.  Finally,  note  how  close  the  x values  are  for  all  the  H atoms, 
falling  in  a small  range  between  0.083  - 0.111.  The  elements  of  the 
H-atom  APTs  are  quite  different  from  each  other,  reflecting  the  differ- 
ent influence  of  the  asymmetric  attached  groups.  A more  complete  and  in 
depth  comparison  could  be  carried  out  to  analyze  the  APTs  according  to 
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contributions  from  charge,  charge-flux  and  overlap  (66).  However,  the 
values  of  the  effective  charge  x give  some  idea  of  how  similar  the 
behavior  of  one  atom  may  be  in  different  molecules. 

Several  sets  of  APTs  for  the  atoms  in  1-Me  U are  obtained  depend- 
ing on  how  we  impose  the  null  condition  (Eq.  (3-5))  to  constrain  to 
zero  the  sum  of  the  APTs.  We  transferred  directly  APTs  for  all  atoms 
of  uracil  (Table  3-3)  except  those  for  and  H y atoms.  For  the  N^-CHg 
fragment  we  transferred  APTs  for  each  one  of  the  sets  in  Table  6-1, 
satisfying  the  constraint  from  the  null  condition  by  choosing  either 
the  APT  for  or  that  for  Cy.  In  the  former  case  we  transferred  the 

APT  for  Cy  listed  in  Table  6-1;  in  the  latter  we  transferred  the  APT 

for  the  atom  from  uracil.  Therefore,  six  complete  sets  of  APTs 
for  1-Me  U were  obtained  following  this  procedure.  Since  the  APTs 
for  the  atoms  in  the  uracil  ring  system  are  given  in  Table  3-3,  only 
the  additional  calculated  APTs  for  the  N-CH2  group  in  1-Me  U are  shown 
in  Table  6-2. 

Let  us  describe  briefly  each  one  of  those  sets  in  Table  6-2: 
set  A uses  the  APTs  from  atom  in  uracil,  the  APTs  of  H atoms  of  the 
N-CH^  group  from  N-methylacetamide,  and  the  APT  for  Cy  atom  was  ob- 
tained by  applying  the  null  condition;  set  B uses  the  APTs  transferred 
from  the  four  atoms  in  the  CH^  group  of  N-CH^  from  N-methylacetamide, 

and  the  APT  for  was  obtained  from  the  null  condition;  set  C uses 

the  APTs  transferred  from  in  uracil  and  from  the  C-CH^  group  in 
N-methylacetamide  for  the  H atoms,  obtaining  the  APT  of  C from  the 
null  condition;  set  D makes  use  of  the  same  APTs  for  the  CH^  atoms 
from  the  C-CH3  group  of  N-methylacetamide,  but  the  APT  of  N^  is  ob- 
tained from  the  null  condition;  set  E makes  use  of  the  APTs  from  N^ 


TABLE  6-2 

SETS  OF  APTs  FOR  1 -METHYLURAC I L (in  e units) 
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dfhe  am  stands  for  N— methylacetamide,  and  am(N— CH3)  and  am(C  CH3)  imply  that  these  APTs  are 
those  of  the  N-CH3  and  C-CH3  groups  of  N-methylacetamide. 

eThe  acetone  APT  was  taken  from  Ref.  39. 

£Px(H15)  is  obtained  by  reversing  the  signs  of  the  off-diagonal  components  of  PX(H14>  involving  z 


112 


in  uracil  and  from  the  H atoms  in  CH3  from  acetone,  with  the  APT  of  C? 
obtained  from  the  null  condition;  set  F uses  all  the  APTs  from  the  CH^ 
group  of  acetone  with  the  null  condition  determining  the  APT  for  the 
atom. 

It  should  be  noted  that  the  pair  of  sets  A and  B,  C and  D,  and 
E and  F are  the  same  except  for  the  null  condition  assumption. 

Therefore,  the  chemical  effect  of  replacing  by  CH^  is  assumed  in 
sets  A,  C,  and  E to  be  concentrated  on  N^;  in  B,  D,  and  F it  is 
assumed  to  be  on  C^.  The  consequences  of  these  assumptions  will  be 
seen  when  we  discuss  the  predicted  intensities  for  the  vibrational 
nodes  of  these  methylated  uracil  molecules. 

3-methyluracil 

The  procedure  for  obtaining  the  APTs  for  3-Me  U to  be  used  in 
Calculation  II  follows  the  one  just  described  for  1-Me  U and  l-Me-3-D  U 
For  the  3-Me  U molecule.  Fig.  6-3  shows  the  atom  numbering  and  geometri 
cal  parameters  used  as  well  as  the  Cartesian  coordinates  for  the  CH^ 
group  in  the  Cartesian  system  of  uracil  (Fig.  3-1).  The  coordinate 
system  for  the  APTs  for  the  atoms  from  Table  6-1  was  then  rotated  by 
+239.30°  (Eq.  (3-4))  to  fit  the  orientation  of  the  N-CH^  group  in 
3-Me  U,  as  shown  in  Fig.  6-4.  The  values  of  the  APTs  in  this  rotated 
coordinate  system  for  3-Me  U are  shown  in  Table  6-3.  Again,  we  ob- 
tained six  sets  of  APTs  for  the  N— CH^  fragment  in  3— Me  U depending  on 
how  the  null  condition  is  satisfied.  These  sets  are  given  in  Table 
6-4.  Their  description  is  analogous  to  that  given  for  the  correspon- 
ding six  sets  of  APTs  for  the  N-CH3  fragment  of  1-Me  U. 
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x 

C?  -2.398326 

H 3 -2.734069 

H, , -2.797948 

14 

H15  -2.797948 


y 

0.090527 

1.127531 

•0.408282 

■0.408282 


z 

0.000000 

0.000000 

0.882944 

-0.882944 


Figure  6-2.  Geometrical  parameters  of  3-Me  U and  Cartesian  coordinates 
for  the  CH3  group  in  the  Cartesian  coordinate  system  of 
uracil.  Coordinates  for  the  ring  atoms  of  uracil  are  given 
in  Table  3-1.  All  C-H  bond  lenghts  in  the  CH3  group  are 
taken  as  1.092  £. 


114 


x 


3 

<U 

Pu 

5S 

3 

1 

O 

r-H 

u 

00 

u 

o 

CO 

44 

3 

CJ 

CO 

1 

H 

Z 

3 

<C 

3 

-C 

0) 

4J 

3 

4-1 

44 

O 

4-1 

O 

CO 

H 

CO 

p- 

CD 

< 

X 

cc) 

CU 

3 

44 

4-1 

o 

• 

3 

3 

e 

•H 

o 

CO 

<0 

•H 

4-1 

X 

4-1 

-O 

i 

CO 

o 

cn 

44 

o 

o 

4-4 

(X 

4J 

o 

-a- 

i 

vO 

0) 

u 

3 

so 

•H 

Pn 


115 


TABLE  6-3 

TRANSFERRED  APTs  (in  e)  FROM  THE  CHg  GROUPS  OF  ACETONE, 
AND  FROM  THE  N-CH3  AND  C-CH3  GROUPS  OF  N-METHYLACETAMIDE , 
AND  THAT  FOR  THE  N3  ATOM  OF  URACIL,  USED  FOR  THE  N-CH3 
GROUP  OF  3-METHYLURACIL 


CH3  of  acetonea,b  N-CH3  of  N-Me  Aa’c  C-CH3  of  N-Me  Aa’c 


-.300 

0.073 

0 1 

0.622 

0.041 

0 ^ 

-.276 

-.134 

0 

w 

0.035 

-.043 

0 

0.199 

0.009 

0 

-.100 

0.176 

0 

0 

0 

0.039, 

0 

0 

0.382 

l 0 

0 

-.010; 

X 

= 0.183 

X 

= 0.437 

X 

= 0.212 

0.042 

0.087 

0 

0.090 

0.077 

0 

0.036 

0.099 

0 

Fx«13> 

0.044 

-.070 

0 

-.011 

0.017 

0 

-.005 

-.143 

0 

0 

0 

0.101 

\ 0 

0 

0.126, 

\ 0 

0 

0.073, 

X 

= 0.094 

X 

= 0.101 

X 

= 0.111 

0.054 

-.041 

0.072 

-.056 

-.072 

0.085^ 

0.049 

-.028 

0.062’ 

Px(H14) 

-.029 

0.040 

0.096 

-.072 

0.048 

0.110 

0.008 

0.072 

0.067 

^0.011 

0.093 

-.028 

0.028 

0.083 

-.080 

.0.019 

0.060 

0.012, 

X 

= 0.102 

X 

= 0.128 

X 

= 0.083 

n3. 

uracil 

n3,  n- 

-Me  A 

c3,  N- 

-Me  A 

'-.623 

-.054 

° \ 

/ -1.016 

-.390 

0 

1.207 

0.385 

0 \ 

^X 

(n3,c3) 

-.094 

-1.403 

0 

-.704 

-.921 

0 

0.356 

1.786 

0 

0 

0 

-.660 

\ 0 

0 

-.815  1 

\ 0 

0 

0.511 

X 

= 0.967 

X 

= 1.032 

X 

= 1.314 

aAll  APTs  have  been  rotated  to  coordinate  axes  of  Fig.  6-2. 
b,c,<^See  footnotes  b,c  and  d of  Table  6-1. 


TABLE  6-4 

SETS  OF  APTs  FOR  3-METHYLURACIL 
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'Notes  b,  c,  d,  and  e of  Table  6-2  apply  to  this  table. 

To  obtain  the  PX(H15)  just  reverse  the  PX(H14)  elements  involving 
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1.3- dimethyluracil 

The  case  of  1,3-di  He  U involves  a combination  of  the  APTs  for 
1-Me  U and  3-Me  U.  First,  Fig.  6-5  shows  the  geometry  of  1,3-di  Me  U, 
where  the  only  thing  we  have  done  is  adding  to  the  1-Me  U geometry 
(Fig.  6-1)  the  CH3  group  geometrical  parameters  taken  from  3-Me  U 
(Fig.  6-3). 

In  order  to  make  the  transfer  of  APTs  in  the  methylated  uracils 
a systematic  and  consistent  procedure,  we  formed  twelve  sets  of  APTs 
for  1,3-di  Me  U combining  sets  A and  B,  C and  D,  and  E and  F for 
1-Me  U (Table  6-2)  with  sets  G and  H,  I and  J,  and  K and  L for  3-Me  U 
(Table  6-4),  respectively.  The  resulting  sets  of  APTs  are  shown  in 
Table  6-5.  This  procedure  guarantees  that  the  sum  of  APTs  closes  to 
zero,  and  also  that  we  are  using  the  same  APTs  for  these  molecules. 
Therefore,  we  can  test  the  transferability  of  APTs  in  this  series 
of  related  compounds.  Of  course,  we  are  expecting  that  the  experi- 
mental and  calculated  intensities,  as  well  as  frequencies,  of 

1.3- di  Me  U to  show  similarities  with  those  of  other  methylated 
uracils . 

Transfer  of  ab  initio  Force  Constants 
In  order  to  carry  out  Calculations  I and  II  we  used  a combination 
of  the  "corrected  to  4-31G"  ab  initio  force  constants  for  uracil  (Ta- 
ble 4-2)  with  force  constants  obtained  for  N-methylformamide , using  a 
4-3 1G  basis  set,  by  Sugawara,  Hirakawa,  Tsuboi,  Kato,  and  Morokuma 
(40).  As  stated  earlier,  the  use  of  ab  initio  force  constants  from  a 
4-31G  basis  set  calculation  normally  reproduces  the  experimental 
frequencies  for  a molecule  after  uniform  scaling. 
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Figure  6-5.  Geometry  and  atom  numbering  for  1 ,3-dimethyluracil 
in  the  Cartesian  coordinate  system  of  uracil. 
Coordinates  for  the  ring  atoms  of  uracil  are  given 
in  Table  3-1,  for  the  N1-CH3  methyl  group  in  Fig. 
6-1,  and  for  the  N3-CH3  methyl  group  in  Fig.  6-2. 


TABLE  6-5 

SETS  OF  APTs  FOR  1 , 3-DIMETHYLURACIL 
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TABLE  6-5  continued 
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The  internal  and  symmetry  coordinates  for  all  four  methylated 
uracil  molecules  considered  in  this  study  are  those  listed  for  uracil 
in  Fig.  4-1  and  Table  4-1,  respectively,  for  Calculation  I (12  atoms). 
1-Me  thy lur ac il 

The  replacement  of  H7  atom  position  in  uracil  (Fig.  3-1)  by  CH3 
for  Calculation  II  adds  9 more  internal  coordinates  to  1-Me  U,  as 
shown  in  Fig.  6-6.  The  new  set  of  in-plane  symmetry  coordinates  for 
1-Me  U is  displayed  in  Table  6-6. 

The  modification  of  the  uracil  force  field  (Table  4-2)  to  include 

the  force  constants  transferred  from  N-methylf ormamide  to  1 Me  U has 

been  underlined  in  Table  6—7.  For  example,  the  diagonal  force  constant 

°-l 

for  the  symmetry  coordinate  (N^-H^  stretch)  of  8.256  md  A in 
uracil  was  changed  to  5.960  md  A ^ (since  it  now  corresponds  to  the 
N^-C^  stretch) . Note  that  Calculation  I is  made  using  only  the  upper 
21  x 21  submatrix  of  Table  6-7,  while  Calculation  II  uses  the  complete 
26  x 26  matrix.  Again,  the  same  force  constants  were  used  for  l-Me-3- 
D U in  Calculations  I and  II. 

3-Methyluracil 

Let  us  now  consider  the  transfer  of  ab  initio  force  constants  for 
3-Me  U.  The  sets  of  in-plane  internal  and  symmetry  coordinates  for 
this  molecule  are  given  in  Fig.  6-7  and  Table  6-8,  respectively.  The 
force  field  for  this  molecule  is  shown  in  Table  6—9,  where  the  modifi- 
cations of  the  uracil  force  field  have  again  been  underlined 
1,3-Dimethyluracil 

For  1,3-di  Me  U,  Fig.  6-8  and  Table  6-10  show  the  in-plane  inter- 
nal and  symmetry  coordinates,  respectively.  The  force  field  which  we 
used  for  this  molecule  was  a combination  of  that  obtained  for  1— Me  U 
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Figure  6-6.  In-plane  internal  coordinates  for  1-Me  U. 
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TABLE  6-6 

IN-PLANE  SYMMETRY  COORDINATES 
FOR  1-Me  U 


S1  -Arl 

vNjC2  stretching 

S2  -Ar2 

vC2N3  stretching 

S3  -4r3 

vNjC^  stretching 

S4  " Ar4 

vC^Cj  stretching 

S5  Ar5 

vCjCg  stretching 

S6  ” Ar6 

vCgNj  stretching 

S?  - 4r? 

vN^c  1 stretching 

S8  ” Ar8 

vC2Og  stretching 

Sg  -*r9 

vN^39  stretching 

S10  " Ar10 

vC^O^q  stretching 

S11  “ Arll 

vCjHjj  stretching 

S12  * Ar12 

vC6H12  strecctlin8 

S13  “ (A813-48u+4B15-4816+4817-4818)//6 

SRing  deformation  I 

S14  ’ < A613-A814-A615+  A816-A617-A618)//n 

SRing  deformation  II 

S15  ■ <A814-A815+A817-A818)/2 

dRing  deformation  III 

S16  - (A819-A825)//5 

dNjC  bending 

S17  * <A820-A826)/,/2 

SCjO  bending 

S18  ’ (A82rA827)/''2 

5N3c  bending 

S19  " (A622-4628)//2 

dC^O  bending 

S20  " (i823~A829^'/,2 

5CjH  bending 

S21  * (A824-A830)/^ 

5CgH  bending 

S22  “ (Ar3i+Ar32+Ar33)/  3 

vMe  stretching 

S23  ” (Ac^+Ac^j+Ac^g-Ao^-Ac^g-Ao^g) / 8 

$Me  symm.  def . 
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Note:  The  symmetry  coordinates  are  defined  in  Table  6-6. 
The  units  are  md  for  stretch-stretch,  md  X for 
bend-bend,  and  md  for  bend-stretch  elements. 
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Figure  6-7,  In-plane  internal  coordinates  for  3-Me  U. 
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TABLE  6-8 

IN-PLANE  SYMMETRY  COORDINATES 
FOR  3 -Me  U 


1 “Arl 

vN^Cj  stretching 

2 *ir2 

vC2N3  stretching 

'3  ’ Ar3 

vN,C,  stretching 
3 4 

'4  “Ar4 

vC.C.  stretching 
4 5 

*5  Ar5 

vC,C,  stretching 
5 b 

*6  " Ar6 

vC.N.  stretching 
6 1 

*7  Ar7 

vN^c  7 stretching 

38  Ar8 

vC208  stretching 

S9  Ar9 

vN^3  9 stretching 

sio  ” Ario 

\iC4010  stretching 

sn  “ Arn 

vC5Hu  stretching 

S12  ” ir12 

vC6H12  stretchin8 

S13  * (A813-A814+AB15-4B16+ABi7-46i8)/'/6 

6Ring  deformation 

I 

S14  ’ ( A813“ABU-A615+  A016-AS17-AS18)//U 

6Ring  deformation 

II 

S15  * (ABu-A815+4617-4S18)/2 

{Ring  deformation 

III 

S16  ’ (AB19-4S25)//5 

6NjC  bending 

3 17  “ (A820-A826)//5 

SCjO  bending 

S18  " (A621-A827)/>/2 

6NjC  bending 

S19  ” (AB22-AS2g)//2 

6C,0  bending 
4 

S20  ” (A823-AB29)//2 

(5C5H  bending 

S21  " (A824_A830)/-^ 

5C,H  bending 
0 

S22  " (Ar3l+Ar32+Ar33) / 3 

vMe  stretching 

S23  - (Aa34+Aa35+Aa36-Aa37-Aa38-Aa39)/  6 

{Me  symm.  def . 

S24  - (2ar31_or32_ctr33^  6 

5Me  deg.  def. 

S25  - (2Aa36“Aa34“Aa35)/  6 

dMe  deg . def . 

S28  - (2Aa37“Aa38~Aa39)  / 6 

rMe  rocking 
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FORCE  CONSTANTS  FOR  3-Me 
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Fig.  6-8.  In-plane  internal  coordinates  for  1,3-di  Me  U. 
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TABLE  6-10 

IN-PLANE  SYMMETRY  COORDINATES 
FOR  1,3-di  Me  U 


T = Ari 

vN^C2  stretching 

;2  = Ar2 

vC2Ng  stretching 

;3  *Ar3 

vN^C^  stretching 

'4  ' Ar4 

vC^Cj  stretching 

'5  Ar5 

vCjCg  stretching 

;6  ■ Ar6 

vCgN^  stretching 

'7  " Ar7 

vMjC  7 stretching 

'8  4r8 

vCjOg  stretching 

'9  “ Ar9 

vN^3g  stretching 

'10  " Ar10 

stretching 

'll  " Arll 

vCjH^  stretching 

;12  * Ar12 

vCgH12  stretching 

>13  ’ (AB13-A814+A615“A616+A617-AS18)/^ 

SRing  deformation  I 

*14  - ( AS13-A614-Aa!5+  A616-AS17-AB18>//ri 

SRing  deformation  II 

'l5  " ^AB14~AB15+AB17-4S18^ /2 

6Ring  deformation  III 

;16  ’ (^19-4625)//5 

<5NjC  bending 

'17  “ (i620-AB26)/,/5 

SCjO  bending 

'is  ■ 

SN^c  bending 

’l9  " C4822-4628)/^ 

SC^O  bending 

'20  “ (a823-462g)/»^ 

6CjH  bending 

'21  " (AB24-AB3C))/^ 

6CgH  bending 

>22  * (Ar3i+Ar32+Ar33)//3 

vMe  stretching 

>23  * (Aa34+Aa35+Aa3g-Aa37-Aa3g-Aa3g) //6 

5Me  symm.  def . 

>24  “ (2ar3i_ar32-ar33)/^ 

5Me  deg.  def. 

>25  “ ( 2Aa36_Aa34~Aa3 5)  //& 

5Me  deg . def . 

>26  “ (2Aa37_Aa38“Aa39)/V6 

rMe  rocking 

>27  “ (Ar40+Ar41+ar42)/^3 

VMe  symm.  def. 

>28  " (A«43+A*44+^5-A*46-A<*47-A<><.48)/7S 

6Me  symm . 

>29  * (2Ar^Q-Ar43-Ar^2) /^*> 

vMe  deg . 

>30  " (2Aai5-Aa44-Aa43)  //6 

SMe  deg. 

>31  = (2Aa46-Aa47-Aoi48)//6 

rMe  rocking 

TABLE  6-11 

FORCE  CONSTANTS  FOR  1,3-di  Me 
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(Table  6-7)  and  that  for  3-Me  U (Table  6-9).  This  force  field  is  given 
in  Table  6-11.  Note  also  that  Calculation  I makes  use  of  only  the  first 
21  x 21  portion  of  this  matrix,  while  Calculation  II  uses  the  complete 
31  x 31  matrix. 

Experimental  and  Predicted  Spectra 

For  the  matrix-isolated  methylated  uracil  molecules  (1-Me  U, 
l_Me-3-D  U,  3-Me  U,  and  1,3-di  Me  U)  the  experimental  infrared  spectra 
have  been  reported  by  Szczesniak,  Nowak,  Szczepaniak,  Chin,  Scott,  and 
Person  (19).  They  interpreted  such  spectra  making  use  of  quantum  me- 
chanical predictions  of  frequencies  and  intensities  of  the  normal  vi- 
brations, with  a calculation  that  is  basically  the  same  as  Calculation  I 
in  the  present  study,  except  that  they  used  the  APTs  for  uracil  reported 
in  Ref.  17,  and  here  we  are  using  the  corrected  APTs  for  uracil  (Table 
3-3).  One  interesting  finding  in  the  study  by  Szczesniak  et  al.  (19) 
is  the  marked  difference  between  spectral  patterns,  both  observed  and 
predicted  in  the  carbonyl  region  for  uracil  and  the  methylated  uracils. 

The  observed  frequencies  and  relative  intensities  of  the  absorption 
bands  given  by  Szczesniak  et  al.  (19)  have  been  collected  and  compared 
with  our  calculated  frequencies  (scaled  by  0.90)  and  absolute  intensi- 
ties obtained  with  Calculations  I and  II,  for  all  the  in-plane  normal 
modes  of  the  methylated  uracils,  as  shown  in  Tables  6-12  to  6-15. 

The  PEDs  are  also  given  in  these  tables  as  part  of  the  assignments  for 
the  in-plane  fundamentals. 

In  order  to  compare  the  experimental  matrix  frequencies  and  inten- 
sities with  the  predicted  results,  we  carried  out  simulations  with  the 
predicted  results  using  Calculations  I and  II,  to  be  compared  with  a 
similar  simulation  of  the  experimental  spectra.  The  simulated  spectra 
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6C.H(29),  6C,H(27),  rock.  Me(25)  1272  1144 
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TABLE  6-13 

RESULTS  FOR  1-METHYL-3-D-URACIL  IN  PLANE-VIBRATIONS 
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TABLE  6-13  continued 
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SNjMe(69),  oC^OC^),  6 Ring  III  (9)  352  316  4 below  studied 

region 

Predicted  Intensity  Sun  (km  mol  *)  : 1952 
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(42),  6 Sing  1(16),  6C_H(14),  vC  N (10)  1351  1216 
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‘All  footnotes  to  Table  6-12  (a  to  1)  apply  to  this  table. 
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RESULTS  FOR  3-METHYLURACIL  IN-PLANE  VIBRATIONS 
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SC6H(26),  511^(22),  vC4C5(14),  <5C5H(13),  vN^CH), 
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Predicted  Intensity  Sum  (km  mol  ) : 1999 
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RESULTS  FOR  1 , 3-DIMETHYLURACIL  IN-PLANE  VIBRATIONS 
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(10),  vC  N (10)  , 744  669 
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for  all  four  methylated  uracils  are  given  in  Figs.  6-9  to  6-12.  All 
simulations  were  done  as  described  in  Chapter  4 for  the  simulated 
spectra  of  uracil. 

We  have  also  collected  in  Tables  6-12  to  6-15  all  the  predicted 
intensities  for  the  methylated  uracils  using  the  different  sets  of 
APTs  given  in  Tables  6— 2,  6—4,  and  6—5.  A close  look  to  these  tables 
shows  that  it  is  difficult  to  decide  which  particular  set  of  APTs 
better  reproduces  the  experimental  intensities.  Some  specific  regions 
of  the  experimental  spectra  are  very  well  reproduced  with  one  set  of 
APTs,  while  the  same  set  of  APTs  does  not  reproduce  other  regions. 

In  order  to  facilitate  the  discussion  of  results,  the  spectra 
have  been  divided  into  three  major  regions:  the  N— H and  C-H  stretching 
region  (3600  - 2000  cm"1),  the  C=0  and  C=C  stretching  region  (2000  - 
1500  cm-'*'),  and  the  ring  stretching  and  deformation  region  (1500  - 
300  cm-*-).  Let  us  now  consider  in  detail  each  one  of  these  spectral 
regions . 

Spectral  Region  Between  3600  - 2000  cm 

From  Tables  6-12  to  6-15  and  Figs.  6-8  to  6-15  we  can  see  that 
for  the  N^H  and  NjH  stretching  vibrations  of  1-Me  U and  3-Me  U the 
observed  frequencies  in  argon  are  somewhat  lower  than  the  predicted 
values,  as  was  seen  for  uracil  (Table  4-3).  However,  the  observed 
wavenumber  difference  between  these  two  N-H  stretches,  about  35  cm  , 
is  close  to  the  predicted  value  of  26  cm  , in  contrast  with  uracil 
where  the  observed  difference  value  was  found  to  be  around  50  cm 
compared  with  the  predicted  value  of  26  cm  1 . The  experimental 
measured  values  (16)  for  the  intensity  of  the  N^H  stretching 
mode  in  3— Me  U vapor  (75  km  mol  ) and  that  of  the  N^H  stretching 
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Figure  6-9.  Simulated  infrared  spectra  for  1-Me  U.  a)  Predicted  spectrum 
using  Calculation  I;  b)  Predicted  spectrum  using  Calculation 
II  with  the  set  of  APTs  F for  acetone  given  in  Table  6-2; 
c)  Experimental  infrared  spectrum  in  an  argon  matrix  (see 
Ref.  19).  The  out-of-plane  modes  are  indicated  by 
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Figure  6-9  continued 
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Figure  6-10.  Simulated  infrared  spectra  for  l-Me-3-D  U.  a)  Predicted 
spectrum  using  Calculation  I;  b)  Predicted  spectrum 
using  Calculation  II  with  the  set  of  APIs  F for  acetone 
given  in  Table  6-2;  c)  Experimental  spectrum  in  an  argon 
matrix  (see  Ref.  19).  The  out-of-plane  modes  are 
indicated  by 
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Figure  6-10  continued 
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Figure  6-11.  Simulated  infrared  spectra  for  3-Me  U.  a)  Predicted 
using  Calculation  I;  b)  Predicted  using  Calculation 
II  with  the  set  of  APTs  L for  acetone  given  in  Table 
6-4;  c)  Experimental  infrared  spectrum  in  an  argon 
matrix  (see  Ref.  19).  The  out-of-plane  modes  are 
indicated  by  . 
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Figure  6-11  continued 
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Figure  6-12.  Simulated  infrared  spectra  for  1,3-di  Me  U. 

a)  Predicted  spectrum  using  Calculation  I; 

b)  Predicted  spectrum  using  Calculation  II 
with  the  set  of  APTs  12  given  in  Table 
6-5;  c)  Experimental  infrared  spectrum 
in  an  argon  matrix  (see  Ref.  19). 
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Figure  6-12  continued 
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mode  of  1-Me  U vapor  (53  km  mol”1)  are  lower  than  the  predicted  values 
of  107  and  98  km  mol-1,  respectively,  but  the  discrepancy  is  still 
within  the  expected  limit  of  a "factor  of  two".  The  experimental 
intensity  ratio  of  the  NH  stretch  in  3-Me  U to  that  of  the  N^H  stretch 
in  1-Me  U is  1.3,  which  is  in  good  agreement  with  the  value  of  1.1  in 
this  study.  The  predicted  value  of  0.7  in  Ref.  17  for  this  ratio  was 
also  due  to  the  use  of  the  incorrect  APTs  for  uracil. 

It  is  also  shown  in  Tables  6-12  to  6-15  and  Figs.  6-9  to  6-12 
that  the  C-H  stretching  modes  of  the  uracil  ring  are  always  predicted 
with  an  intensity  of  about  10  km  mol  1,  or  10  % that  for  N-H  stretches, 
as  it  was  for  the  uracil  molecule  (Table  4-3).  Experimentally,  only 
a trace  of  absorbance  is  observed  for  such  C-H  stretches  in  the  matrix- 
isolated  spectra.  However,  for  the  C-H  stretches  of  the  methyl  groups 
on  the  methylated  uracils.  Calculation  II  predicts  different  intensities, 
which  depend  on  the  set  of  APTs  used  (Tables  6-12  to  6-15) . The  lowest 
predicted  intensity  values  are  for  the  case  when  the  APTs  of  the  H atoms 
of  acetone  are  used,  which  is  in  better  agreement  with  experimental 
data,  because  these  absorptions  are  not  observed  in  the  matrix-isolated 
spectra.  The  APTs  of  the  H atoms  of  N-methylacetamide  do  not  predict 
good  intensity  values  for  these  C-H  stretches. 

Spectral  Region  Between  2000  - 1500  cm 

The  experimental  spectra  for  the  methylated  uracils  show  a 
complex  structure  of  the  bands,  because  more  than  the  expected  bands 
are  observed  in  the  matrix  (19).  An  inspection  of  Tables  6-12  to 
6-15  and  Figs.  6-9  to  6-17  shows  that  the  frequencies  of  carbonyl 
(C20  and  C40)  bands  of  methylated  uracils  are  predicted  to  be  slightly 
higher,  and  the  separation  between  them  predicted  lower,  than  the 
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observed  frequencies  and  separations  of  the  bands.  It  will  be 
interesting  to  compare  the  experimental  and  predicted  sums  and  ratios 
for  the  intensities  of  modes  6 (C20  stretching)  and  7 (C^O  stretching) 
using  Calculations  I and  II.  This  comparison  is  shown  in  Table  6-16 
for  1-Me  U,  l-Me-3-D  U,  3-Me  U,  and  for  1,3-di  Me  U.  For  1-Me  U 
the  calculated  ratio  of  the  intensity  of  mode  6 to  that  of  mode  7 
in  Calculation  I (APTs  of  uracil)  is  3.47,  which  is  about  1.8  times 
the  experimental  value  of  1.93.  Also,  we  can  see  from  Table  6-16  that 
the  different  sets  of  APTs  used  in  Calculation  II  give,  in  general, 
lower  ratios,  but  still  larger  than  the  experimental  value.  In 
uracil  (Table  4-3)  the  ratio  of  the  intensity  for  both  stretches  is 
2.34  (experimentally)  and  2.17  (predicted).  The  decrease  of  the 
experimental  ratio  going  from  uracil  (2.34),  to  1-Me  U (1.93), 
to  l-Me-3-D  U (1.25),  to  3-Me  U (0.34),  to  1,3-di  Me  U (0.48)  has 
been  discussed  by  Szczesniak,  Nowak,  Szczepaniak,  Chin,  Scott,  and 
Person  (19) . 

For  mode  8 (C=C  stretching)  the  observed  and  scaled  frequency 
values  agree  quite  well  in  all  the  methylated  uracils;  however,  the 
intensity  prediction  for  this  mode  is  always,  using  either  Calcula- 
tions I or  II,  much  lower  than  the  experimental  values. 

Spectral  Region  Between  1500  - 300  cm  ^ 

An  examination  of  Tables  6-12  to  6-15  and  Figs.  6-9  to  6-12 
shows  that  there  are  only  a few  absorption  bands  which  have  similar 
frequencies  and  intensities  in  the  experimental  spectra  of  methyla- 
ted uracils.  For  example,  one  of  these  is  the  band  near  1490  cm  , 
which  according  to  our  PEDs  calculation  has  predominant  contribution 
from  the  CH^  degenerate  deformation  vibration.  In  order  to  interpret 


TABLE  6-16 

COMPARISON  OF  RESULTS  FOR  MODES  v&  AND  v?  OF  THE  METHYLATED  URACILS 
mode  Expl.  Calcn.  I Calcn.  II 
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this  region,  it  is  very  important  to  carry  out  Calculation  II,  since 
here  is  where  we  find  all  the  in-plane  deformation  modes  of  the  methyl 
groups.  Moreover,  because  in  Calculation  II  we  are  using  different 
sets  of  APIs  we  can  compare  the  predicted  intensities  obtained  mainly 
because  methyl  group  deformations.  In  the  particular  case  of  monomeric 
1-Me  U the  assignment  for  most  of  the  bands,  based  on  our  calculations 
given  in  Table  6-12  and  Fig.  6-9,  agrees  quite  well  with  previous 
assignments  reported  in  the  literature  (1,4,13).  However,  the  corre- 
spondence is  not  entirely  straightforward  because  the  previous  spectra 
(1,4,13)  are  related  to  hydrogen  bonded  molecules  in  the  crystal, 
while  our  spectrum  corresponds  to  the  monomeric  molecules,  only 
slightly  perturbed  by  the  interaction  with  the  matrix.  A close  loox 
to  Table  6-12  and  Fig.  6-9  shows  that  the  results  of  Calculations  I 
and  II  do  not  reproduce  the  experimental  spectrum  in  this  region  as 
we  expected,  in  particular  the  "rocking"  deformation  mode.  The  same 
pattern  is  observed  for  the  other  methylated  uracil  molecules  (Tables 
6-13  to  6-15  and  Figs.  6-10  to  6-12).  It  is  also  interesting  to  note 
that  in  the  methyl  deformation  region  the  experimental  intensity  is 
better  reproduced  when  the  APTs  of  the  N— CH^  group  of  N— methylacetamide 
are  used. 

The  1358  cm“l  band  in  the  spectrum  of  1-Me  U is  not  present  in  the 
spectra  of  the  other  methylated  uracils.  It  is  predicted  at  1394  cm 
as  having  an  intensity  of  10  km  mol  ^ (Calculation  I)  and  1393  cm 
with  an  intensity  ranging  from  10  to  26  km  mol  ^ (Calculation  II) . 

It  is  assigned  to  the  NH^  bending  mode. 
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Summary 

We  have  seen  in  this  Chapter  that  vibrational  parameters  (force 
constants  and  atomic  polar  tensors)  obtained  at  the  ab  initio  level 
with  a 4-3 1G  basis  set  for  uracil,  acetone,  N-methylf ormamide,  and 
N-methylacetamide  can  be  transferred  to  predict  infrared  frequency  and 
intensity  data  of  the  methylated  uracil  molecules.  These  ab  initio 
quantum  mechanical  calculations  using  a 4-31G  basis  set  predict  within 
a f actor-of-two  the  experimental  absolute  intensities  and  within  +20 
cm-^  the  experimental  frequencies.  These  calculations,  even  in  this 
limited  form,  are  quantitative  enough  to  be  very  useful  to  assign  and 
to  interpret  the  complex  infrared  spectra  of  these  molecular  systems. 
Furthermore,  we  can  have  considerable  confidence  in  similar  predictions 
for  similar  molecules. 

We  should  certainly  be  able  to  use  the  ab  initio  force  constants 
and  APTs  presented  in  this  dissertation  to  predict  intensities  accurate 
to  a factor-of-two  for  other  pyrimidines,  like  cytosine.  Therefore, 
the  technique  of  transferring  vibrational  parameters  to  predict  the 
infrared  spectra  of  other  molecules  is  a very  promising  one. 

Since  the  calculated  results  presented  in  this  dissertation  are 
just  a first  step,  future  research  in  this  theoretical  and  experimental 
area  of  biomolecuies  could  be  directed  to  the  development  of  more 
refined  schemes  for  transferring  molecular  parameters,  with  the  aim 
of  predicting  more  accurate  intensity  and  frequency  data  which 
reproduces  much  more  closer  the  experimental  data.  Perhaps,  the 
next  step  to  follow  with  this  type  of  calculation  could  be  the 
application  of  a perturbation  theory  refinement  procedure  to 
achieve  this  objective.  Of  course,  more  experimental  work  should 
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be  done  employing  matrix-isolation  techniques,  including:  infrared. 


Fourier  transform  infrared,  Raman  and  resonance  Raman,  to  get  more 
data  of  biological  molecules. 
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